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1.   INTRODUCTION 

Vapor  bubble  nucleation  and  growth  phenomena  have  aCCracted  the 
attention  of  researchers  for  many  years.   Nearly  all  of  the  attention 
has  been  placed  on  saturated  nucleate  boiling  and  the  transition  to  film 
boiling.   Studies  of  boiling  under  transient  conditions  have  been 
devoted  primarily  to  investigations  of  the  transition  to  film  boiling. 
The  nature  of  boiling  phenomena  preceding  the  transition  has  received 
little  attention.   However,  these  phenomena  are  important  aspects  of  the 
broader  problem  of  understanding  the  sequence  of  events  leading  to 
boiling  crises  in  high-pressure  heat  exchangers  such  as  nuclear  power 
reactors. 

Transient  initiation  of  boiling  has  been  described  by  steady-state 
nucleation  models,  which  use  the  size  of  active  nucleation  sites  to 
determine  superheats  necessary  for  bubble  nucleation  from  these  sites. 
The  required  superheats  for  a  simple  transient  situation  can  be  achieved 
by  an  increase  in  the  heater  surface  temperature  or  by  a  decrease  in  the 
fluid  saturation  temperature  by  means  of  a  reduction  in  the  fluid 
pressure. 

Surface  characteristics  are  important  in  the  description  of  the 
sizes  of  cavities  available  as  nucleation  sites.   The  pressure  and 
temperature  conditions  experienced  prior  to  the  initiation  of  boiling 
were  also  postulated  to  be  of  major  importance  by  Fabic  (1),  since 
potential  sites  could  be  flooded  by  pre-pressurization.   Studies 
performed  to  examine  transient  boiling  caused  by  decompression  (2-5) 
have  employed  conditions  such  that  bulk  liquid  flashing  occurs  by  the 
end  of  the  decompression.   Consequently,  the  behavior  when  boiling 


initiation  conditions  are  not  met  for  the  predicted  cavity  sizes  and  the 
reactivation  of  flooded  surface  cavities  has  not  been  thoroughly 
investigated. 

In  this  work,  the  transition  from  convection  to  boiling  during 
pressure  transients  was  investigated  for  a  horizontal,  cylindrical 
platinum  heating  element  submersed  in  a  fluid.   Maximum  overpressures 
were  applied  to  the  system  before  subjecting  the  heating  element  to 
constant  power  delivery,  maintained  during  transients  by  a  custom  made 
power  supply  driven  by  12  volt  wet  cells.   Separate  runs  were  necessary 
to  measure  the  system  pressure  as  a  function  of  time  and  the  element 
temperature  as  a  function  of  time.   The  element  temperature  was 
determined  using  resistance  thermometry.   The  data  were  recorded  using  a 
digital  oscilloscope,  stored  on  a  floppy  disk,  then  transferred  to  a 
computer  for  analysis.   The  time  at  which  the  first  bubble  was  seen  or 
heard  was  recorded  as  the  measured  boiling  initiation  time.   Initiation 
times  were  also  estimated  based  on  the  point  where  the  analyzed  heater 
temperature  began  to  decrease.   Although  these  two  methods  did  not 
always  produce  equivalent  time  values,  ranges  for  true  boiling 
initiation  times  were  determined  from  these  measurements. 

Pressure  transients  displayed  exponential  reduction  behavior  over 
most  of  the  pressure  decrease,  and  they  were  reproducible  in  that 
region.   Boiling  initiation  times  provided  conditions  which  tended  to 
support  the  cavity  deactivation  by  pre-pressurization  hypothesis  of 
Fabic,  although  the  effects  of  initial  pressures  were  not  as  pronounced 
as  expected  by  the  model.   Boiling  was  observed  even  in  cases  where  none 
was  predicted,  but  subsequent  nucleation  site  reactivation  and  recovery 
of  the  heater  temperature  to  steady-state  was  impeded  by  increases  in 
maximum  pressure  applied. 


2.   REVIEW  OF  PREVIOaS  WORK 

2.1  Convection  Heat  Transfer 

In  a  system  composed  of  a  body  immersed  in  a  fluid,  convection  heat 
transfer  will  occur  if  the  body  and  the  bulk  fluid  temperature  are  at 
different  temperatures.   Quantifying  the  heat  transfer  analytically  is 
not  a  simple  matter  because  of  the  dependence  not  only  on  numerous  fluid 
properties,  but  also  on  flow  conditions  and  surface  geometry. 

Boundary-layer  theory  has  been  used  to  derive  numerical  solutions 
for  horizontal  cylinders.   These  solutions  have  been  shown  to  be  quite 
accurate  for  moderate  Rayleigh  numbers  characterising  cases  where  the 
wake  is  confined  to  a  small  region  at  the  rear  of  the  cylinder. 
Churchill  and  Chu  (6)  examined  a  wide  range  of  data  on  steady-state 
convective  heat  transfer  from  horizontal  cylinders.  They  recommended 
the  following  correlation  as  being  applicable  to  both  constant  cylinder 
surface  temperature  and  constant  heat  flux  for  all  Rayleigh  and  Prandtl 
numbers : 

Nu^  =  0.60  +  0.387  f(Pr)Ra^^^  (2.1) 

where  the  Rayleigh  and  Nusselt  numbers  are  based  on  cylinder  diameter, 
and 

f(Pr)  =  [1  +  (0.559/Pr)^^^^]"^''^^  (2.2) 

Morgan  (7)  also  carried  out  a  thorough  analysis  of  published 
experimental  data  on  steady-state  convection  from  horizontal  cylinders, 
and  recommended  the  following  correlation  as  being  applicable  for  gases 
and  liquids  (except  liquid  metals) : 

(2.3) 


where  values  of  B  and  m  are  given  in  the  following  table: 


Ra 

B 

m 

10-1°  _  jo-2 

0.675 

0.058 

10-2  _  j„.2 

1.02 

0.148 

10^  -  10^ 

0.850 

0.188 

10^  -  10^ 

0.480 

0.250 

10^  -  ioi2 

0.125 

0.333 

Fujii,  et  al.  (8)  reported  that  for  convection  from  horizontal 
cylinders,  the  following  correlation  shows  good  agreement  with  both 
calculated  and  measured  heat  transfer  rates: 

2/Nu  =  Jin  [1  +  4.918g(Pr)Ra-"],  (2.4) 


where 


and 


g(Pr)  =  [1  +  (0.492/Pr)^^l^]  *^'  (2.5) 


n  =  0.25  +  1.0/(10  +  5Ra°'l^^).  (2.6) 


Table  2.1  shows  a  comparison  of  these  three  correlations  for  the 
Raylelgh  number  range  of  interest  in  this  work. 

2.2  Boiling  Initiation 

2.2.1  Steady-State  Boiling  Initiation 

Many  theoretical  analyses  and  interpretations  of  experimental  data 
on  boiling  initiation  have  been  based  on  concepts  introduced  by  Griffith 
and  Wallis  (9)  and  embellished  by  Hsu  (10),  Bergles  and  Rohsenow  (11), 


Table  2,1.   Comparison  of  predictions  of  Churchill-Chu, 
Fujii,  and  Morgan  correlations  for  natural 
convection  heat  transfer  from  a  0.25  mm 
diametric  cylinder.   Ambient  conditions: 
T  =  20°C.   Fluid  properties  are  evaluated 
at  the  mean  temperature  between  T  and  T  , 


C) 

Heat 

Flux  (kW/m 

=  ) 

T  -T  (° 
w  a 

Churchi 

.11- 

-Chu 

Morgan 

Fujii,  et  al. 

10 

26. 

2 

29;  4 

31.1 

20 

60. 

3 

68.0 

72.5 

30 

100 

113 

121 

40 

145 

163 

176 

50 

195 

218 

237 

60 

249 

278 

304 

70 

308 

341 

376 

80 

371 

409 

453 

90 

438 

483 

535 

100 

509 

564 

622 

Han  and  Griffith  (12),  Madejski  (13),  and  Cole  and  co-workers  (14,  15, 
16).   In  its  simplest  form,  the  concept  may  be  described  as  follows. 

Consider  present  within  a  fluid  a  heating  element  within  the 
surface  of  which  is  an  active  boiling  initiation  site  with  a  circular 
aperture  of  radius  r   (Fig.  2,1).   As  the  temperature  of  the  heating 
element  is  increased  beyond  the  saturation  temperature  of  the  fluid,  a 
vapor  bubble  emerges  from  the  aperture.   The  surface  of  the  bubble  is 
spherical  shaped  due  to  the  effect  of  surface  tension.   The  minimum 

radius  of  curvature  of  the  bubble  is  r  ,  the  radius  of  the  aperture.   At 

c  *^ 

this  minimum  radius,  surface  forces  are  maximum,  and  the  excess  of  the 
vapor  pressure  within  the  bubble  over  the  static  fluid  pressure  must  be 
maximum. 

If  the  fluid  surrounding  the  bubble  is  at  temperature  T  ,  the  excess 
pressure  is  given  by 

In  which  T   Is  the  saturation  temperature  of  the  fluid  at  system  pressure 
p  =  p  (T  ) ,  p  (T  )  is  the  vapor  pressure  of  the  fluid  at  temperature  T  , 
and  a  is  the  fluid-vapor  surface  tension  evaluated  at  temperature  T  . 
The  excess  temperature  of  the  fluid  at  the  fluid-vapor  interface  of  the 
bubble  can  be  related  to  the  excess  pressure  in  the  bubble  through  use 
of  the  Clauslus-Clapeyron  relation.   Application  of  this  relation  leads 
to  the  approximation  that 


2aT 

f,(r  )  E  T  +  -; — ?-  .  (2.8) 

1   c     s   Xp  r 


where  p   is  the  vapor  phase  density,  assumed  to  be  much  less  than  that  of 
the  liquid  phase,  and  X  is  the  latent  heat  of  vaporization.   This  is  the 


j^    ■■•     ^-^r 


'■>r 


FIG.  2.1.   Illustration  of  bubble  growth  from  an  active 
nucleation  site. 
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threshold  temperature  which  must  be  exceeded  for  the  bubble  to  grow 

beyond  radius  r  . 

Now  suppose  that  the  fluid  surrounding  the  heating  element  has  a 

temperature  distribution  T(r)  =  f  (r) ,  a  function  of  distance  r  from  the 

surface.   Let  T  and  T  represent  the  surface  temperature  and  the  ambient 
w      a   '^  ^ 

temperature  respectively.   If  q  is  the  heat  flux  from  the  surface  to  the 
fluid,  then  approximately 

f^(r)  =  dT/dr  =  -q/k,  (2.9) 

and 

f„(r)  =  T(r)  =  T  -  qr/k.  (2.10) 

i.  w 

where  k  is  the  thermal  conductivity  of  the  fluid. 

Griffith  and  Wallis  argued  that  for  a  bubble  to  grow  in  size  beyond 

the  radius  r  ,  the  surface  temperature  of  the  heating  element  must 

exceed  the  threshold  temperature,  i.e.,  that  f^(0)  =  f,(r  ).   Hsu 

contended  that  for  the  bubble  to  grow,  it  is  necessary  that  the  fluid 

temperature,  not  at  the  heating-element  surface  or  base  of  the  bubble, 

but  at  the  cap  of  the  bubble,  exceed  the  threshold,  i.e.,  that  f„(r  )  = 

z  c 

f  (r  ).   Han  and  Griffith  argued  that  neither  the  temperature  at  the 

bubble  base  nor  that  at  the  bubble  cap  controlled  the  bubble  growth. 

Rather,  it  was  reasoned,  the  average  liquid  temperature  surrounding  the 

bubble  was  the  controlling  factor.   This  led  to  the  criterion  f^(ar  )  = 

f . (r  ),  where  it  was  argued  that  a  =  3/2.   Sakurai  and  Shiotsu  (17) 

accommodated  in  their  arguments  an  empirical  relationship  between  the 

cavity  radius  and  the  bubble  radius  at  boiling  Inception,  namely  that 

the  bubble  radius  is  br  .   Allowing  for  this  emperlcisim,  the  boiling 

initiation  criterion  may  be  stated  as  f^Car  )  =  f , (br  ),  or 

z   c     1   c 


1 
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qar    2aT 

V   c 

For  the  Griffith  and  Wallls  approach,  a  =  0  and  b  =  1.   For  that  of  Hsu, 
a  =  1.6  and  b  =  1.25.   For  that  of  Han  and  Griffith,  a  =  1.5  and  b  =  1. 
Sakarai  and  Shlotsu  recommended  that  a  =  1.8  and  b  =  1.57. 

Suppose  now  that  there  are  present  on  the  heating  element  surface 
nucleatlon  sites  of  a  wide  variety  of  radii.   Cole  (18),  following  the 
prescription  of  Bergles  and  Rohsenow  (11),  showed  that,  as  the  heat  flux 
or  surface  temperature  is  Increased,  the  first  nucleatlon  site  to  become 
active,  i.e.,  from  which  the  first  bubble  will  emerge,  is  that  one  for 
which  r^  =  r^  determined  by  f'(ar^)  =  f;(br^).   From  Eqs .  (2.8)  and 
(2.9), 

r^  =  (2kaT^/p^Xb  =  q)'^.  (2.12) 

Substitution  of  this  value  into  Eq.  (2.11)  leads  to  the  following 
expression  for  the  minimum  surface  superheat  required  for  boiling  from  a 
surface  with  a  full  size  range  of  nucleatlon  sites: 

(T  -  T  )  .   =  (8aT  q/kXp  )'^.  (2.13) 

w    s  min       s^   "^v 

Bergles  and  Rohsenow  offered  the  following  correlation  of  computed  values 
for  water  over  the  pressure  range  15  to  2000  psla: 

q  =  1.55  X  10"*  p^-^"  (1.8(T^-T^)]2-0^7''P°'°^^''    (2.14) 

where  q  is  measured  W/m^,  T  is  K,  and  p  is  MPa.   If,  on  the  other  hand, 

the  largest  nucleatlon  site  present  has  r  =  r   <  r.,  then  the  minimum 

cm* 

surface  superheat  required  for  boiling  initiation  is  the  greater  value 


10 


qar    2aT 

(T  -T  )  ,  =  —-B.  +  —^—  ,  (2.15) 

w  s  min    k     Xp  br 


Selected  values  of  minimum  boiling  superheat  for  H„0  are  given  in  Table 

2.2.   For  each  saturation  temperature,  the  last  four  values  of  r   listed 

m 

are  the  r^^  values  for  the  heat  fluxes  listed.  Thus,  superheats  above  and 

to  the  left  of  the  horizontal  lines  are  for  the  cases  r   <  r.. 

ra    * 

One  may  ask  how  the  superheat  at  boiling  inception  is  affected  by 
the  subcooling  of  the  ambient  liquid.   Suppose  that  a  full  range  of 
initiation  sites  is  available,  so  that  the  boiling  superheat  is  governed 
by  Eq.  (2.13).   Suppose  that  natural  convection  is  the  heat  transfer 
mechanism,  so  that 


q  =  C(T^-T^)"  (2.16) 

where  n  ^   1.2.   By  equating  the  heat  fluxes  in  Eqs.  (2.13)  and 

(2.16),  and  noting  that  T-T  =T  -T  +AT  ^,  one  can  show  that  (T  - 
w  a    w    s     sub  w 

T  )  .   is  given  by  solution  of  the  equation 
s  min    o  J  1 

(T  -  T  )  .   =  C'[(T  -  T  )  .   +  AT  ^l"''^        (2.17) 
w    s  mm        w    s  min     sub 

I, 
where  C  =  (8oT  C/kXp  )^.   Evaluation  of  (T  -T  )  .   is  illustrated  in 
s     V  w  s  min 

Fig.  2.2  where  it  may  be  seen  that,  as  AT  ,  increases,  so  also  does  the 

sub 

minimum  boiling  superheat.   For  example,  for  water  at  T  =  100°C,  an 

increase  in  AT  ,  from  10  to  30°C  causes  an  increase  of  about  S'C  in  the 
sub 

minimum  boiling  superheat. 

Madejski  (13)  developed  a  model  for  boiling  initiation  which  allows 
for  a  non-spherical  bubble  shape  and  avoids  specification  of  the  distance 
from  the  heating  surface  at  which  the  superheat  criterion  must  be 
applied.   As  corrected  by  Schmidt  and  Cole  (14),  the  result  of  the 
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Table  2.2.   Minimum  boiling  superheat  as  a  function  of  maximum  cavity 
radius. 


T  (°C)      r   (urn)      q(MW/m^) 


100        1 
2 

4.73 
6.69 

10.6 
15.0 

150        0.5 
1 

2.31 
3.26 
S.16 
7.30 

195       0.5 
1 

1.31 
1.86 
2.93 
4.15 


0.1 

0.2 

0.5 

1.0 

32.9 

33.1 

33.5 

34.3 

16.7 
7.6 
5.9 
4.7 
4.4 

17.0 
8.3 
6.9 
6.2 
6.2 

17.9 

10.4 

9.8 

9.8 

9.8 

19.3 
13.8 
13.8 
13.8 
13.8 

15.4 

15.5 

15.7 

16.1 

7.8 
3.7 
2.8 

2.2 

8.0 
4.0 
3.3 

3.0 

8.4 
5.0 
4.7 

4.7 

9.1 
6.6 
6.6 

6.6 

2.1 

3.0 

4.7 

6.6 

5.1 

5.2 

5.4 

5.8 

2.7 
2.1 

1.6 

1.3 

2.8 
2.3 

1.9 

1.7 

3.3 
2.9 

2.7 
2.7 

4.0 
3.9 
3.9 
3.9 

1.2 

1.7 

2.7 

3.9 

12 


log  (q) 


EQ.  (2.16) 


Eq.  (2.17),  ^T^ub>  0 


log'VTs) 


FIG.  2.2.   Determination  of  boiling  Initiation  superheat 
requirements  based  on  Eqs.  (2.13)  and  (2.16). 
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analysis  is  expressed  as  the  ratio  of  the  minimum  surface  superheat  in  a 
variable  temperature  situation  to  that  required  if  the  liquid  were 
uniformly  superheated ,  namely 


w"  s  min  ,  ^^^^^^  ^^^  _^  ^/^^^  ^2.18) 


(T  -T  )      """"^  m"'  w 

W   S  w 


The  function  increases  monotonically  with  its  argument,  reaching  a  value 
of  about  4.5  at  an  argument  of  3.   The  utility  of  the  Madejski  approach 
is  limited  by  the  availability  of  data  for  boiling  initiation 
temperatures  in  uniformly  superheated  liquids, 

2.2.2  Transient  Boiling  Initiation 

Most  of  the  work  performed  concerning  transient  boiling  phenomena 
has  focused  on  constant  pressure  systems  subjected  to  transient  heating. 
Few  are  of  interest  here  because  pressure  effects  generally  were  not 
considered.   Those  studies  that  used  pressure  as  a  parameter  (e.g., 
references  19,  20,  and  21)  reported  that  the  superheat  required  for 
boiling  initiation  decreased  with  increasing  pressure.   However,  these 
studies  examined  boiling  at  various  constant  system  pressures  and  did  not 
test  the  effects  of  applying  a  temporary  overpressure  before  an 
experimental  run. 

Fabic  (1)  thoroughly  examined  the  consequences  of  the 

pressure-temperature  history  of  a  heating  element  on  boiling  initiation. 

He  argued  that  the  maximum  radius  of  an  active  (unflooded)  initiation 

site  at  the  time  of  a  test  is  governed  by  the  maximum  value  of  the  ratio 

(p~P  )/cJ  experienced  prior  to  the  test.   Here  p  is  the  system  pressure 

and  p  the  fluid  vapor  pressure.   Call  this  maximum  ratio  (p  -p   ) /a  . 
V  '^o  vo   o 

The  maximum  cavity  radius  is 
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cos(ir-9  ) 

r  =  2cr  —  (2.19) 

c     o  p  -p 


where  0   (>90°)  is  the  contact  angle  between  the  cavity  wall  and  the 
vapor-liquid  interface  within  the  cavity.   The  value  of  9   is  taken  as 
108"  for  water  at  all  temperatures  and  pressures.   Values  of  the  ratio 
and  r  are  given  in  Table  2.3  for  water  under  various  conditions. 

Fabic  noted  that  most  test  fluids  at  some  time  are  at  20°C  and 
atmospheric  pressure.   This  leads  to  an  upper  bound  of  r  =  r  =0.5  ym, 
for  radii  of  nucleation  cavities.   He  argued  that  boiling  initiation 
takes  place  with  a  contact  angle  9  between  0  and  66",  and  at  a  surface 
temperature  T  which  satisfies  the  equation 

■  p  (T  )  -  p  =  2acos9  /r* ,        t*<t    ,  (2.20) 

v  w  r  c     c-  c 

=  2acos9  /r  ,    r'^>r  ,  (2.21) 

r  c     c   c 

where  a   is  evaluated  at  T  .   Upper  and  lower  limits  on  T  would  depend  on 

the  choice  of  0°  or  66°  for  9  .   If,  in  fact,  the  largest  potential 

initiation  site  had  a  radius  r  <  r  ,  then  that  radius  would  apply. 

m    c  '^'^  ' 

One  may  inquire  of  the  effect  of  boiling  from  the  heating  surface  at 

test  conditions  for  purposes  of  degassing  the  surface.   Fabic  argued  that 

the  process  would  "arm"  previously  inactive  (flooded)  sites  with  vapor 

and  thus  that  the  pressure  and  temperature  under  which  the  degassing 

took  place  would  define  r  . 
c 

One  may  also  inquire  of  the  effect  of  degassing  at  saturation 
conditions.   Fabic  is  silent  on  this  point.   Presumably  sites  of  all 
sizes  would  be  activated  and  r  would  govern  subsequent  boiling 
initiation. 
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Table  2,3,   Maximum  active  site  radius  as  a  function  of 
temperature  and  pressure,  computed  using  the 
Fabic  Model  [1], 

Fluid  Temp.  (°C)    Pressure  (MPa)     (p-p  ) /o  (um"^)     r'^  (ym) 
V c 

20  0.101 
1.398 

70  0.101 

90  0.101 

95  0.101 

100  0.476 

100  1.123 

100  1.398 

120  0.476 

140  0.476 

145  0.476 

165  1.398 

185  1.398 

190  1.398 


1.35 

0.46 

19.1 

0.032 

1.08 

0,57 

0.51 

1.21 

0.28 

2.22 

6.33 

0.098 

17.3 

0,036 

21.9 

0.028 

5.02 

0.12 

2.24 

0.28 

1.21 

0.51 

15.2 

0.041 

6.64 

0.093 

3.55 

0.17 
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Winterton  (2)  also  developed  an  expression  for  conditions  at 
nucleation  based  on  the  pressure-temperature  history  of  the  heater, 

p  +  p  -  p  = (p  -  P   -  P   )  »  (2.22) 

V   *^g   ^  y  a      ^^o   Wo   '^go 


where  p   is  the  partial  pressure  of  noncondensible  gas,  and  y  is  defined 


cos  (tt  -  0  +  -t) 

Y  = —   ,  (2.23) 

cos(e^  -  j) 

where  a  is  the  included  apex  angle  of  the  cavity.   It  will  be  noted  that 
Eq.  (2.22)  reduces  to  Eq.  (2.20)  for  the  case  of  a  well  degassed  system 
and  a  cylindrical  cavity  (a  =  0) ,   Winterton  performed  a  number  of 
experiments  using  vertical  tubes  to  test  theories  of  nucleation  by 
surface  cavities.   He  examined  both  transient  heat  flux  -  constant 
pressure  and  constant  heat  flux  -  transient  pressure  cases.   Liquid 
pressures  were  in  the  range  of  approximately  0.01  to  0.1  MPa  for  both 
cases  and  true  heat  fluxes  in  the  transient  heat  tests  were  10  kW/m^ 
maximum,  resulting  in  a  rate  of  rise  in  temperature  of  less  than  l°C/s. 
For  the  transient  heat  case,  he  reported  that  the  vapor  pressure  at 
nucleation  increased  with  increasing  overpressures,  but  that  the  effect 
was  not  as  pronounced  as  indicated  by  the  theory  with  y  set  equal  to 
unity.   He  reported  that  the  order  of  magnitude  of  y  was  1,  but  his 
experiments  indicated  that  the  value  was  not  constant.   Instead,  y 
increased  as  the  overpressure  increased.   The  limitations  of  his 
experimental  apparatus  precluded  making  conclusions  about  the  effect  on 
the  transient  pressure  case.   However,  for  the  conditions  under  which  he 
was  able  to  make  both  transient  heat  and  transient  pressure  measurements, 
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he  found  no  significant  difference  in  the  calculated  cavity  radius 
required  for  nucleation,  thus  confirming  the  equivalence  of  the  two 
methods  of  nucleation,   Gallagher  and  Winterton  (4)  performed  a  study 
which  expanded  on  the  power  transient  boiling  work  of  Winterton.   Again, 
a  very  low  heat  flux  was  used,  and  overpressures  up  to  0.7  MPa  were 
applied  before  performing  tests  at  0,0122  MPa.   They  found  that  special 
preparation  of  the  test  section  was  necessary  to  produce  consistent 
results.   They  reasoned  that  this  preparation  was  associated  with 
obtaining  suitable  and  consistant  values  of  contact  angles  within  the 
surface  cavities.   Their  results  show  a  strong  effect  of  pressure  on 
deactivating  cavities.   They  reported  good  agreement  between  their 
results  and  Eq,  (2.22)  using  a  constant  value  for  Y  of  5,4. 

Faw  and  VanVleet  (22),  investigating  boiling  initiation  from  a 
horizontal  wire  in  a  subcooled  pool  of  water,  reported  initiation 
conditions  that  were  not  consistent  with  any  single  size-distribution  of 
nucleation  sites.   However,  nearly  all  results  were  consistent  with 
predictions  made  using  Fabic's  model.   They  also  found  that 
pre-pressurization  well  above  saturation  pressure  led  to  greatly 
increased  superheat  requirements  at  boiling  initiation  consistent  with 
the  model.   They  noted  that  only  the  first  test  after  pressurization  was 
so  affected;  subsequent  tests  in  the  series  had  boiling  initiation 
superheats  approximately  those  expected  in  the  absence  of 
pre-pressurization.   They  postulated  that  the  first  test  reactivated 
previously  flooded  sites. 

Investigation  of  boiling  phenomena  during  transient  pressure 
conditions  has  been  prompted  by  desire  to  determine  the  response  of  a 
high  pressure,  high  temperature  heat  exchanger  system  (e.g.,  a  nuclear 
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reactor)  to  a  loss  of  pressure.   Since  the  mechanical  integrity  of  such  a 
system  can  be  breached  by  events  resulting  from  inadequate  heat  transfer, 
most  studies  have  examined  the  transition  from  nucleate  to  film  boiling 
and  critical  heat  fluxes  (e.g.,  references  23,  24,  and  25). 

The  initiation  of  boiling  during  decompression  has  received  limited 
attention.   Fable  argued  that  his  pressure-temperature  history  modeling 
of  boiling  initiation  would  also  apply  to  the  pressure  transient  case. 
He  noted  that  the  nucleation  event  could  occur  due  to  a  decrease  in 
system  pressure  even  If  the  surface  temperature  T  did  not  change.   In 
this  case  the  system  pressure  p  necessary  for  vapor  nucleation  could  be 
predicted  from  Eqs.  (2.20)  and  (2.21).   Values  of  the  boiling  Initiation 
pressure  in  water,  based  on  the  Fable  model,  are  given  in  Table  2.4  for 
water  under  various  conditions. 

It  can  be  expected  from  Fable's  model  that  an  Increase  in  the 
initial  overpressure  applied  to  a  system  should  cause  an  increase  in  the 
Instantaneous  surface  superheat  [T  -T  (p) ]  necessary  to  initiate 
boiling.   This  Increase  would  manifest  Itself  by  producing  a  measurable 
delay  between  the  time  boiling  initiation  is  predicted  by  steady-state 
considerations  and  when  it  is  actually  observed.   Hooper  and  Abdelmessih 
(26)  observed  just  such  an  occurrence  in  a  study  of  liquid  flashing. 
However,  Kenning  and  Thirunavukkarasu  (27)  found  no  indication  of  delay 
in  bubble  growth  in  their  study  of  bubble  nucleation  characteristics  of 
surfaces. 

Weisman,  et  al.  (3),  in  a  study  of  boiling  initiation  using  heated 
and  unheated  ribbons,  found  that  measurable  delays  in  boiling  initiation 
could  be  encountered  under  appropriate  conditions.   They  suggested  that 
the  results  of  Kenning  and  Thirunavukkarasu  could  be  explained  by  the 
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Table  2.4.   Boiling  initiation  system 

pressure  values  prescribed  by 
the  Fabic  model  (1)  at  T  = 
lOCC  and  using  6  =  66°f 
T  =  T  until  after  initial 
pressure  is  applied. 


P,.  (MPa)  T  CO  p   (MPa) 

1  WD 


0.239  140            0.204 

0.273  0.165 

0.308  0.125 

0.308  150            0.250 

0.377  0.175 

0.411  0.137 

0.377  160            0.330 

0.446  0.258 

0.515  0.186 

0.584  0.U4 

0.515  170            0.380 

0.584  0.312 

0.653  0.243 

0.790  0.106 

0.584  180            0.546 

0.653  0.481 

0.790  0.350 

0.997  0.155 


'  1 
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fact  that  the  pressure  was  released  very  rapidly  to  well  below  saturation 
pressure,  thus  activating  unflooded  cavities  and  allowing  bubble  growth 
to  begin  almost  immediately.   In  their  own  investigation,  the  authors 
also  found  that  the  superheat  required  to  initiate  boiling  on  an  unheated 
ribbon  increased  with  increasing  rate  of  depressurization.   This  is 
consistent  with  results  obtained  in  studies  of  the  response  of  hot  water 
to  rapid  depressurization  (28) (29) .   Furthermore,  Weisman  et  al.  found 
that  at  low  decompression  rates  the  superheats  required  to  initiate 
boiling  showed  no  effects  due  to  the  pressure  history,  and  instead 
approached  steady-state  values.   These  superheats,  in  the  3  -  4°C  range, 
were  substantially  below  those  predicted  for  the  pressure-time  histories 
considered.   The  authors  reasoned  that  small  bubbles  of  noncondensible 
gas  could  remain  at  a  cavity  base  after  application  of  an  overpressure. 
As  pressure  is  reduced,  vapor  diffuses  into  the  bubble,  which  then  fills 
the  cavity.   Thus,  an  active  cavity  is  present  with  a  radius  larger  than 
the  maximum  unflooded  radius  predicted,  and  boiling  is  possible  at  a 
superheat  below  that  predicted  using  the  maximum  unflooded  radius.   The 
authors  speculated  that  the  finite  time  to  fill  such  a  cavity  could  be 
the  reason  for  the  delay.   Heated  ribbons  were  used  to  verify  that  the 
first  boiling  in  the  system  occurred  on  the  ribbon.   The  authors  noted 
that  the  superheats  were  higher  than  for  the  tests  without  heating,  but 
they  did  not  comment  on  how  these  superheats  compared  to  predicted 
values. 

Sakurai,  et  al.  (5),  examined  transient  boiling  of  a  test  heater 
caused  by  depressurization  with  near-exponential  reduction  periods 
ranging  from  3  to  60  ms.   The  platinum  wire  test  heater  was  initially  in 
a  nonboiling  state  in  a  pool  of  water,  and  the  heat  generation  rate  at 
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the  heater  was  held  constant  throughout  the  experimental  run.   They  used 

initial  pressures  of  0.59,  1.08,  and  1.9  MPa  with  water  temperatures  of 

353  and  373  K.   Rapid  decompression  was  achieved  using  a  rupture  disc 

device  to  vent  the  test  section  to  atmosphere.   The  authors  found  that 

the  heater  surface  temperature  and  true  heat  flux  remained  constant  until 

boiling  initiated.   The  incipient  boiling  point  was  assumed  to  have  been 

reached  at  the  time  the  heater  temperature  began  to  decrease.   They 

argued  that  this  gave  the  criterion  for  boiling  initiation  even  if  it  may 

not  have  been  the  true  incipient  boiling  point.   They  found  that  the 

pressure  p  ,  and  the  instantaneous  surface  superheat,  AT,  ,  at  the 
D  b 

incipient  boiling  point  depended  upon  the  heat  flux,  even  if  the  initial 
overpressure  and  the  reduction  period  were  the  same.   They  noted  that  AT 
increased  with  decreasing  p,  and  was  little  dependent  on  the  pressure 
reduction  period,  consistent  with  the  incipient  boiling  model  of  Sakurai 
and  Shiotsu  [Eq.  (2.11)]. 

2.3  Nucleate  Boiling 

A  well-known  correlation  for  fully-developed  subcooled  and  saturated 
nucleate  pool  boiling  is  that  of  Rohsenow  (30) ,  which  may  be  expressed  in 
the  form 

q  =  Ci(T„  -  T^)^  (2.24) 

The  factor  C-  is  a  function  of  pressure  (fluid  properties  evaluated 
under  saturation  conditions)  and  fluid-surface  combination,  namely 

C^  =  [(PfP^Jg/o]*^  yXPr^"  (c/XC^)^,  (2.25) 

where  g  is  the  acceleration  of  gravity,  u  is  the  viscosity  of  the  fluid, 
and  c  is  the  heat  capacity  of  the  fluid.   The  constant  n  and  the  value 
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of  C„  depend  on  the  fluid-surface  combination  (18),  and  take  the  values 
1.0  and  0.013  respectively  for  a  platinum  surface  in  water.   Unless 

otherwise  indicated,  properties  are  those  of  the  liquid  phase.   For 

-2  -3 
example,  for  platinum-water  at  atmospheric  pressure,  C.  =  140  Wm  K 

Stephan  and  Abdelsalam  (31)  carried  out  an  exhaustive  investigation 

of  experimental  data  on  pool  boiling  in  various  fluids  and  recommended 

the  following  correlation  as  being  applicable  to  water 

q  =  C3  (1^-1^)2-°'^  (2.26) 

in  which  C„  is  a  function  of  saturation  temperature.   Based  on  data  in 
the  paper,  C_  may  be  expressed  empirically,  in  S.I,  units,  as 

C^  =  229  p^'^'^^  (2.27) 

where  p  is  in  units  of  MPa  for  the  range  0.1  to  20.   As  an  example,  for 

—  7  —"^   Oft 
boiling  at  atmospheric  pressure,  C^  =  52  Wm  K   '   .   Thus,  for  a  given 

surface  superheat,  this  correlation  predicts  a  much  lower  heat  flux  than 

that  predicted  by  the  Rohsenow  correlation. 
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3.   THEORY 

3.1  Determination  of  Largest  Potentially  Active  Cavity 

In  the  absence  of  free  gas  bubbles,  nucleation  sites  for  boiling 
initiation  will  be  located  at  surface  cavities  which  have  vapor  trapped 
in  them.   The  rate  of  boiling  heat  transfer  is  governed  by  the  spatial 
distribution  and  the  size  spectrum  of  nucleation  sites.   Heat  flux  and 
superheat  required  to  initiate  boiling  generally  depend  on  the  size  of 
the  largest  site.   The  largest  potential  nucleation  site,  i.e.,  a 
cavity,  idealized  by  a  conical  shape,  containing  vapor,  can  be 
determined  by  considering  the  pressure-temperature  history  of  the 
surface  along  with  the  fluid-surface  wetting  characteristics. 

Consider  an  arbitrary  surface  cavity,  idealized  by  a  conical  shape, 
containing  vapor  with  a  fluid  of  infinite  extent  covering  the  surface 
(Fig.  3.1).  Assume  a  well  degassed  fluid  so  that  the  partial  pressure 
of  noncondensible  gas  is  negligible  compared  to  vapor  pressure  and  fluid 
pressure.  The  system  is  initially  at  pressure  p-.  and  at  temperature  T. 
at  which  the  vapor  pressure  is  p  .  less  than  p_.  and  the  surface  tension 
is  a.  . 

When  a  pressure  p-„  >>  p  .  is  applied,  the  fluid  will  enter  the 
mouth  of  the  cavity.   Mechanical  equilibrium  across  the  meniscus 
requires 

2a 

where  R-,  is  the  radius  of  curvature  of  the  meniscus  (see  Fig.  3.1). 
From  a  geometric  consideration,  the  cavity  radius  r  at  the  meniscus  line 
of  contact  is  given  by 
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FLUID 


FIG.  3.1.   Illustration  of  vapor  entrapment  at 
potential  nucleation  site. 
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r  =  R^  cosCtt  -   e^  +  -), 


2a 

cos(iT   -   e      +  !■)  (3.2) 


Pf2  -  Pvl  ^2 


where  6   is  the  advancing  contact  angle  between  the  fluid  and  the  cavity 

wall,  and  a  is  the  included  angle  of  the  cavity  apex.   The  minimum 

pressure  necessary  for  the  meniscus  to  enter  a  cavity  with  mouth  radius 

r  can  also  be  obtained  from  Eq.  (3.2).   If  this  requirement  is  not  met, 

then  the  contact  angle  is  less  than  6  and  the  meniscus  is  at  the  mouth 

a 

or  the  cavity  is  flooded. 

The  argument  of  the  cosine  must  fall  between  0  and  Sj'"  since  r  must 
exist  for  vapor  to  be  trapped  in  the  cavity.   This  provides  a  maximum 
value  for  the  apex  angle  of  the  cavity, 

a    <  29  -  TT.  (3.3) 

max     a 

If  a  is  greater  than  this  value,  the  cavity  will  flood  with  fluid  and 

become  deactivated.   It  is  also  apparent  from  Eq.  (3.3)  that  6  must  be 

greater  than  ^^tt  for  vapor  to  be  trapped  in  the  cavity. 

Equations  (3.1),  (3.2),  and  (3.3)  describe  the  potential  nucleation 

sites  available  for  the  pressure-temperature  history  of  the  system  when 

p.,  equals  p,    .   Cavities  with  mouth  radii  greater  than  or  equal  to 
1/         t,max  ^ 

the  critical  radius  r  =  r  when  P-^  =  p    and  apex  angles  less  than 
c      '^f2   "^max      "^     " 

a    will  have  the  meniscus  located  at  the  position  inside  the  cavity 
max  '^  ^ 

where  the  local  cavity  radius  is  equal  to  r  .   Cavities  with  apex  angles 

greater  than  a    will  be  flooded,  while  cavities  with  mouth  radii  less 
max 

than  r  will  not  be  affected. 
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The  description  of  the  meniscus  behavior  as  conditions  leading  to 
nucleation  are  imposed  depends  on  the  manner  in  which  these  conditions 
arise.   Boiling  can  be  initiated  by  decreasing  the  fluid  pressure, 
increasing  the  vapor  pressure  by  heating,  or  a  combination  of  these  two 
actions. 

3.2  Activation  of  a  Nucleation  Site 

3.2.1  Decrease  in  Fluid  Pressure 

If  the  fluid  pressure  is  released  to  a  value  p^  such  that  p  ,  <^ 
p   <  pp^,  the  radius  of  the  meniscus  will  change  to  a  new  value  given  by 

The  radius  of  the  meniscus  line  of  contact  remains  constant  at  the  value 
given  by  Eq.  (3.2),  and  the  contact  angle  6  must  satisfy 


2a 

^     cos(7r  -  6  +  ^)  =  r  .  (3.5) 


Pf  -  Pvl  2 


As  P-.  approaches  p  ,  ,  the  contact  angle  9  must  decrease  so  that  r 
remains  constant.   However,  as  long  as  p-  is  greater  than  p  .,  the 
meniscus  will  be  concave  on  the  fluid  side. 

When  p-  becomes  less  than  p  , ,  the  meniscus  radius  becomes  negative 
by  Eq .  (3.4),  i.e.,  the  meniscus  flips  from  concave  to  convex  on  the 
fluid  side.   From  the  force  balance. 


K   = .  (3.6) 
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With  no  boiling,  r  remains  constant  at  the  value  determined  by  Eq, 
(3.2).   However,  r  can  also  be  defined  by 

r  =  ^   cos(8  -  f) 


2a 

cos(e  -  f)  (3.7) 

Pf  -  Pvl         ^ 


Equation  (3.7)  shows  that  r  will  remain  constant  as  6  decreases  to 
the  value  of  the  receding  contact  angle  6  .   When  9  equals  9   the  fluid 
will  begin  to  recede  over  the  previously  wetted  cavity  wall  and  the 
meniscus  line  of  contact  will  move  towards  the  cavity  mouth.   This  will 
occur  when 


Pvi-Pf  =7^"^  (^-f^-  "-^^ 


When  the  meniscus  is  at  the  cavity  mouth,  boiling  inception  will  take 
place  when  the  meniscus  radius  of  curvature  is  at  a  minimum,  i.e.,  at 
the  cavity  mouth  radius,  r  .   This  will  occur  at  a  pressure  p"  such  that 

2a 

p„,  -  py  =T-*-  (3.9) 


The  application  of  these  equations  to  a  range  of  cavity  sizes  will  be 
discussed  later.   However,  it  is  apparent  that  p"  <  pi  when  r  =  r  . 

3.2.2  Increase  in  Vapor  Pressure 

If  a  temperature  T  >  T  is  applied  at  pressure  p.,  where  p  .  < 
Pj  <  Pf2'  ^^^  surface  tension  and  the  vapor  pressure  will  assume  new 
values  a,  and  p  „,  respectively.   If  T   is  less  than  the  saturation 
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temperature  at  p   the  meniscus  radius  of  curvature  will  change,  but  the 
line  of  contact  will  not,  and  the  radius  r  at  the  line  of  contact  will 
remain  constant  at  the  value  given  by  Eq.  (3.2).   Here, 


2°2 
K, .  (3.10) 

^       Pf  -  Pv2 


and 


r  =  Rj^  cos(Tr  -   9   +  y) 


2a, 

cosdr  -   e   +  ^)    =  r    .  (3.11) 


Pf   -  Pv2  2c 


Again,  the  Interface  is  stable  with  the  meniscus  concave  on  the  fluid 
side  for  (p^  -  p  „)  >  0. 

When  the  applied  temperature  T-  becomes  greater  than  the  saturation 
temperature  at  p_,  the  meniscus  radius  of  curvature  becomes  negative  by 
Eq.  (3.10)  and  the  meniscus  flips  from  concave  to  convex  on  the  fluid 
side.   The  force  balance  in  this  case  requires 


2^2 


The  radius  of  the  meniscus  line  of  contact  is  given  by 
r  =  R^  cos(e  -  f) 


2°2 

cos(6  -  ^)  =  r  .  (3.13) 


Pf  -  Pv2         2 


The  meniscus  radius  of  contact  remains  constant  as  6  decreases  to  the 
value  of  the  receding  contact  angle  6  .   Movement  of  the  meniscus  line 
of  contact  towards  the  cavity  mouth  will  commence  when 
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c 


corresponding  to  a  temperature  T',   Again,  boiling  will  initiate  when 

the  meniscus  radius  of  curvature  reaches  a  minimum  equal  to  the  cavity 

mouth  radius  r  .   This  will  occur  at  a  temperature  T"  such  that 
c  z 

2a'' 

P;2-Pf=-^-  ^^-^5) 

r 
c 

Although  a  decreases  with  increasing  temperature,  under  usual  conditions 

p".  >  p'„  or  T''  >  Tl  when  r  =  r„. 
v2    v2     2    2       c    c 

3.3  Cavity  Geometry 

Consider  now  a  broad  size  range  of  active  cavities  present  on  the 

surface.   The  system  has  experienced  conditions  such  that  the  meniscus 

radius  of  curvature  is  convex  on  the  fluid  side.   If  the  meniscus  line 

of  contact  is  below  the  cavity  mouth,  an  increase  in  p  ~  P^r  can  not  be 

balanced  in  Eq.  (3.8)  or  (3.1A),  and  the  line  of  contact  will  move 

towards  the  cavity  mouth.   First  boiling  will  then  occur  from  a  cavity 

with  mouth  radius  r  such  that 
c 

r  =  r  sec(8  -  f)  ,  (3.16) 

c    c      r   2 

since  the  threshold  conditions  for  meniscus  motion  will  also  satisfy  the 
requirement  for  nucleation  from  a  cavity  with  mouth  radius  r  . 

The  effect  of  the  cavity  apex  angle  in  determining  potential 
nucleation  sites  is  revealed  by  Eq.  (3.2),  which  shows  that  vapor  can  be 
trapped  at  lower  liquid  temperatures  as  a  decreases.   In  other  words, 
given  two  cavities  with  the  same  mouth  radius,  the  cavity  with  the 
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smaller  apex  angle  will  boil  first.   Therefore,  a  cylindrical  cavity  (o 
=  0)  will  be  the  largest  potentially  active  nucleation  site  available, 
and  the  superheat  required  to  achieve  the  receding  contact  angle  9 
will  be  less  than  for  a  conical  cavity.   The  meniscus  will  rush  out  of 
the  cylindrical  cavity  for  an  increase  in  p  -  Pf.  and  the  relation  of 
Eq.  (3.9)  or  (3.15)  must  be  satisfied  to  initiate  boiling  if  the 
entrance  angle  to  the  cavity  raouth  is  90°.   The  superheat  in  this  case 
may  be  higher  than  that  required  for  a  conical  cavity  with  mouth  radius 
r  if  a  is  sufficiently  small.   However,  Fable  (1)  argued  that  the 
cavity  mouth  would  be  rounded,  and  appropriate  rounding  at  the  entrance 
would  allow  boiling  to  initiate  after  the  meniscus  flip  with  no 
additional  superheat.   This  special  rounding  would  also  apply  to  conical 
cavities,  but  the  only  effect  would  be  to  allow  the  superheat  necessary 
for  boiling  inception  in  a  cavity  with  mouth  radius  r  to  approach  the 
value  for  r  .   Since  Fable's  argument  provides  the  lower  limit  for 
boiling  initiation,  further  discussion  will  consider  cylindrical 
cavities  where  r  =  r  is  the  radius  of  the  largest  active  cylindrical 
cavity  according  to  Eq.  (3.2). 

3.4  Effect  of  Contact  Angles 

The  contact  angles  used  in  this  work  are  the  same  as  those  adopted 

by  Fabic.   He  selected  an  advancing  contact  angle  6  =  108°  for  water  on 

metal  for  all  values  of  system  pressures  and  temperatures.   This  value 

is  based  on  reported  findings  that  the  largest  contact  angle  water  can 

make  with  a  solid  surface  is  6  =  105°  to  U0°.   His  selection  of  the 

a 

receding  contact  angle  8  =  66°  was  more  arbitrary  since  it  was  a  best 
fit  value  for  his  work.   However,  it  also  agreed  with  the  observance  of 
a  40°  maximum  difference  between  6  and  6  . 
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Consider  a  horizontal  surface  immersed  in  well  degassed  water  of 
infinite  extent.   The  surface  and  the  fluid  are  initially  at  uniform 
temperature  T.  when  an  initial  pressure  p,„  is  applied.   The  surface  is 
then  heated  to  temperature  T  =  T  such  that  no  boiling  occurs. 
Beginning  at  time  t  =  0,  the  system  is  subjected  to  a  time-dependent 
decrease  in  the  fluid  pressure  p-(t) .   No  boiling  will  occur  until  p, 
reaches  the  value  necessary  for  boiling  initiation. 

The  radius  of  the  largest  active  cylindrical  cavity  is  obtained 
from  Eq.  (3.2)  for  the  conditions  T.  and  p., 

2°1 
-'     -  cos(iT  -  e  ).  (3.17) 


Since  r  remains  constant  until  boiling  is  initiated,  Eq.  (3.13)  also 
applies , 


2°2 
cos(e  ).  (3.18) 


^v2 


The  fluid  pressure  necessary  for  boiling  initiation,  p,  ,  can  be  obtained 

b 


by  combining  Eqs.  (3.14)  and  (3.15), 

cos(9  ) 


Pb  =  Pv2  -  ^  cos(7r  -   6  )  (Pf2  "  Pvl^  "  ^^'^^ 


Figure  3.2  shows  the  effect  of  surface  temperature  on  the  incipient 

boiling  pressure,  where  it  is  observed  that  increasing  T   causes  p,  to 

w        ^b 

increase.   The  figure  also  shows  that  for  appropriate  extreme  conditions 
of  initial  pressure,  no  boiling  is  predicted. 
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The  Influence  of  the  advancing  and  receding  contact  angles  on  the 
boiling  initiation  pressure  can  be  examined  using  Eq.  (3.19). 
Increasing  the  ratio  of  the  cosine  terms,  by  decreasing  6  or  increasing 
6  ,  requires  a  lower  value  of  p   (and  hence  a  higher  superheat)  to 
satisfy  the  condition  for  nucleation.   Decreasing  the  ratio  of  the 
cosine  terms  has  the  opposite  effect.   The  effect  of  selected  changes  in 
contact  angles  on  the  ratio  of  the  cosine  terms  is  shown  in  the 
following  table,  and  the  effect  on  p   for  a  representative  case  is  shown 
in  Fig.  3.3.   It  is  apparent  that  the  value  of  p  is  very  sensitive  to 
small  changes  in  the  value  of  either  contact  angle. 


g 

e 

r 

cos (6  ) 

a 

cos(Tr  -  6  ) 

108° 

66° 

1.316 

110° 

63.25° 

1.316 

106° 

68.73° 

1.316 

106° 

66° 

1.476 

108° 

62.86° 

1.476 

108° 

68.45° 

1.189 

110° 

66° 

1.189 
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4.   APPARATUS  AND  EXPERIMENTAL  PROCEDURE 

4,1  Pressure  System 

4.1.1  Pressure  Vessel 

The  test  vessel,  shown  in  Fig.  4.1,  is  made  of  type  316  stainless 
steel.   The  vessel  is  40.64  cm  in  length  with  an  outer  diameter  of  12.065 
cm  and  an  inner  diameter  of  9.8425  cm.   The  end  flanges  are  sealed  with 
Ultek  oxygen  free  hardened  copper  0-rings  Model  2  68-4000,  from 
Perkin-Elmer  Corp.   One  end  flange  is  penetrated  by  a  2  kW  Firerod 
electrical  heater.  Model  N6A22,  from  Watlow,  Inc.   This  heater  is  used  to 
maintain  the  ambient  temperature  of  the  test  fluid.   The  other  end  flange 
is  penetrated  by  the  electrode  assembly.   A  stainless  steel-sheathed 
copper-constantan  thermocouple  from  Watlow,  Inc.  penetrates  the  cylinder 
wall  to  provide  constant  monitoring  of  the  bulk  fluid  temperature. 

The  two  optical  ports  are  essential  to  the  observation  of  bubble 
growth  behavior  on  the  electrically  heated  test  element.   The  design  of 
the  window  housing  is  shown  in  Fig.  4.2.   The  windows  are  of  T08 
commercial  grade  clear  fused  quartz  obtained  from  Amersil,  Inc.   The 
discs  are  4.445  cm  in  diameter  and  1.5  cm  thick  and,  as  estimated  by 
Amersil,  a  pressure  of  13.6  MPa  may  be  used  with  an  approximate  safety 
factor  of  3.5.   The  cushions  and  seals  for  the  windows  are  fabricated  of 
teflon. 

4.1.2  Pressure  Measurement  System 

A  Kistler  Quartz  Miniature  Pressure  Transducer  Model  603A  penetrates 
the  cylinder  wall  to  provide  the  means  to  measure  the  pressure  history 
during  a  transient.   The  point  of  penetration  is  between  one  of  the 
optical  ports  and  the  end  flange  containing  the  electrode  assembly,  and 
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THERMOCOUPLE 


FIG.  4.1.   Schematic  diagram  of  the  pressure  vessel. 
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TEFLON    WASHER 


TEFLON    RING    .6  3  5  cm 
WIDE    WITH    TAPERED    END 


QUARTZ    WINDOW 


FIG.    4.2.      Pressure  vessel  window  housing. 


38 

on  the  horizontal  plane  through  the  test  element.   This  location  places 
the  pressure  transducer  in  closest  proximity  to  the  test  element  so  that 
recorded  information  represents  the  actual  pressure  history  of  the  test 
element  during  the  experiment.   Although  the  transducer  will  withstand 
sustained  temperatures  of  up  to  260°C  without  damage,  its  temperature 
must  be  limited  to  about  120°C  to  minimize  drift  and  ensure  good  low 
frequency  response.   A  water  cooling  adapter  designed  by  Kistler  provides 
this  cooling  capability,  as  well  as  being  a  most  convenient  means  of 
attachment  to  the  pressure  vessel.   Figure  4.3  shows  details  of  the 
cylinder  wall  where  the  adaptor  is  mounted,  and  Fig.  4,4  shows  the 
cooling  adaptor  design  from  Kistler.   Since  this  design  is  for  a 
physically  larger  transducer,  an  adaptor  was  designed  to  mount  the  Model 
603A  pressure  transducer  in  the  cooling  jacket.   The  design  of  the  second 
adaptor  is  shown  in  Fig.  4.5.   Both  adaptors  were  fabricated  by  machine 
shop  personnel  in  the  Department  of  Physics. 

The  quartz  sensing  element  in  the  pressure  transducer  is  oriented 
such  that  increasing  applied  pressure  results  in  a  negative-going  charge 
signal  that  is  inverted  to  a  positive-going  voltage  signal  by  Kistler 
charge  amplifiers.   The  Kistler  Model  504E  Dual  Mode  Amplifier  is 
designed  to  develop  the  full  capabilities  of  the  quartz  transducer  while 
eliminating  the  attenuating  effects  of  the  cable  and  transducer 
capacitance. 

Three  factors  must  be  considered  when  setting  up  the  pressure 
measurement  system  for  optimum  performance:  the  system  transfer  function 
(i,e,,  volts  output  per  measurand  unit  input),  the  correct  time  constant 
for  the  event,  and  the  amplifier  drift.   All  are  directly  affected  by  the 
settings  on  the  dual  mode  amplifier. 
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FIG.  4.3.   Pressure  transducer  cooling  adaptor 
mounting  to  the  pressure  vessel. 
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FIG.  4.4.   Pressure  transducer  cooling  adaptor,  designed  by  Kistler. 
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CONNECTOR  ADAPTOR 
5/  I  6"  HEX 


M7  x  0.7S  THD. 


MODEL  603A  PRESSURE 
TRANSDUCER 


RUBBER  SEAL  RTNG 


MATERIAL-TYPE  304  STAIN 
LESS  STEEL 


FIG.  4.5.   Pressure  transducer  adaptor  for  mounting  the 
Model  503A  to  the  cooling  adaptor. 


42 

The  system  transfer  function  is  inversely  proportional  to  the 
transducer  sensitivity,  which  is  set  on  the  dual  mode  amplifier  as  a 
value  between  1.00  and  11.00.   Since  the  pressure  transducer  used  in  this 
study  has  a  sensitivity  less  than  1.00  pC/psi,  a  range  multiplier  of  10 
is  automatically  required  in  determining  the  transfer  function.   The 
transfer  function  is  also  dependent  on  the  amplifier  gain,  which  varies 
directly  with  the  range  switch  position. 

The  time  constant  is  selected  according  to  the  temporal  behavior  of 
the  input  signal.   Since  the  depressurization  rates  used  in  this  study 
were  small,  the  long  time  constant  setting  was  chosen  to  avoid  signal 
decay  before  the  end  of  the  transient.   Unfortunately,  the  Improvement  in 
low  frequency  response  allowed  by  the  long  time  constant  setting  is  at 
the  expense  of  amplifer  drift  as  a  function  of  time. 

Amplifier  drift  is  inversely  proportional  to  the  range  setting.   For 
the  long  time  constant  the  drift  is  measured  in  millivolts  per  second. 
Thus  the  need  to  minimize  drift  must  be  balanced  against  the  desire  to 
establish  a  system  transfer  function  that  provides  maximum  deflection  of 
the  output  signal  for  a  given  pressure  change. 

4.1.3  Pressurizer 

The  pressurization  system,  as  shown  in  Fig.  4.6,  consists  of  a 
bottle  of  commercial  nitrogen,  the  pressurizer.  Model  30A-10WS,  from 
Greer  Hydraulics,  Inc.  and  the  pressure  vessel.   Seamless  stainless  steel 
tubing  (0.635  cm)  with  0.635-cm  taper-seal  valves  from  High  Pressure 
Equipment  Company  is  used  for  gas  delivery.   The  high  pressure  nitrogen 
enters  one  side  of  the  pressurizer  and  inflates  a  neoprene  rubber  bag. 
This  bag  presses  against  the  water  in  the  other  half  of  the  pressurizer. 
Thus,  pressure  can  be  applied  to  the  test  fluid  without  having  the 
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pressurizing  gas  dissolve  in  the  test  fluid.   The  pressure  vessel  is  also 
connected  to  the  water  side  of  the  pressurizer  using  seamless  stainless 
steel  tubing.   On  the  gas  side,  there  are  two  Asco  Red  Hat  solenoid 
valves.  Model  8262A214  from  Automatic  Switch  Co.,  one  to  control  inflow 
of  gas  and  another  solenoid  valve  to  control  the  flow  of  gas  to  the 
atmosphere.   On  the  fluid  side  of  the  system,  there  is  another  solenoid 
valve.   Once  the  system  is  pressurized,  this  valve  can  be  closed  so  that 
the  pressure  in  the  pressure  vessel  remains  constant  while  the  pressure 
on  the  gas  side  of  the  bladder  is  raised  or  lowered.   Then  the  valve  can 
be  opened  causing  an  instantaneous  compression  or  decompression.   There 
are  also  two  Model  63-5622  pressure  gauges  from  Matheson  Corp.  and  one 
Model  7108P-200  pressure  gauge  from  Transcat  associated  with  the  system. 
One  Matheson  gauge  is  on  the  pressure  vessel  itself  and  the  other  is  on 
the  supply  line  of  the  nitrogen  gas,  while  the  Transcat  gauge  is  on  the 
gas  side  of  the  pressurizer. 

4.2  Heating  Element 

The  test  element,  shown  in  Fig.  4.7,  consists  of  two  brass  rods  each 
0.6350  cm  in  diameter  and  a  teflon  block.   The  teflon  block  is  used  to 
prevent  the  rods  from  sagging.   The  vertical  support  assembly  can  be 
adjusted  to  provide  for  different  lengths  of  platinum  wire.   The  teflon 
seals  at  the  penetration  points  of  the  electrodes  served  a  double  purpose 
of  providing  a  water  tight  seal  and  also  providing  insulation  between  the 
electrodes  and  the  pressure  vessel. 

The  heating  element  is  of  pure  platinum  wire  (SPPL-010  from  Omega 
Engineering,  Inc.)  of  diameter  2.50  x  10~  m  and  is  mounted  horizontally. 
The  maximum  length  of  the  heating  element  is  9.6  cm.   The  spring 
attachment  provides  tension  across  the  wire.   The  platinum  wire  not  only 
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acts  as  the  test  element  but  also  as  a  temperature  sensor.   The 
resistance,  R(6),  of  pure  platinum  has  a  well  characterized  temperature 
dependence,  (32) 

R(6)  =  Rjj  (1  +  ae  +  B9  =  ), 
where 

a  =  3.92  X  10"^ 
6  =  -5.5  X  10"^ 

9  =  the  temperature  of  the  platinum  wire  in°C, 
R(9)  =  the  resistance  of  the  wire  at  temperature  6  in  ohms, 
R.  =  the  resistance  of  the  platinum  wire  at  0°C. 
Once  the  resistance  is  measured,  the  average  temperature  of  the  platinum 
wire  can  be  calculated. 

4.3  Electrical  System 

The  electrical  system  consists  of  two  twelve  volt  storage  batteries, 
a  control  system,  and  a  voltage  bias  system.   The  two  storage  batteries 
are  wired  in  parallel  and  provide  the  dc  current  to  the  test  element. 
The  control  system  regulates  current  to  the  test  element  such  that  the 
power  delivered  remains  constant  with  time.   Figure  4.8  is  a  schematic 
diagram  of  the  control  system  circuit.   The  voltage  bias  system  is  an  EAI 
Analog  Computer,  Model  TR-10.   Simple  linear  summing  circuits  apply 
bucking  voltages  to  the  test  voltages  across  the  standard  resistor  and 
the  test  element  so  that  only  the  variations  in  these  voltages  are 
recorded  by  the  oscilloscope. 

4.4  Data  Measurement  System 

Data  for  the  experimental  runs  are  taken  on  a  Nicolet  Explorer  III 
Digital  Oscilloscope  equipped  with  a  floppy  disc  storage  system  and  a 
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FIG.  4.8.   Control  system  and  power  amplifier  schematic  diagrams. 
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RS-232C  serial  interface  port  for  transmission  of  the  stored  data  to 
other  devices.   The  oscilloscope  has  two  inputs,  each  with  separate 
amplifying  circuits.   The  oscilloscope  is  set  in  Cursor  mode  for 
triggering.   In  this  mode  the  trigger  source  signal  is  continuously 
monitored  until  the  triggering  conditions  are  met  or  exceeded.   At  this 
instant  the  designated  storage  channel  is  assigned  trigger  time  zero,  and 
input  signal  regions  both  preceding  and  following  the  trigger  time  are 
stored  for  observation.   Thus,  although  the  event  time  zero  may  not 
correspond  to  the  trigger  time  zero,  the  entire  transient  can  be  captured 
by  judicious  choice  of  the  storage  channel  for  trigger  time  zero.   Any 
difference  between  trigger  and  event  time  zero  can  be  accounted  for  in 
analysis . 

Three  voltage  signals  are  available  during  a  transient:  one  relating 
to  pressure  vessel  pressure,  and  two  relating  to  test  element 
temperature.   Experimental  runs  to  determine  the  pressure  history  during 
decompression  are  performed  by  recording  the  voltage  signal  from  the 
pressure  transducer  amplifier,  with  the  oscilloscope  triggered  by  this 
signal.   The  unused  input  channel  is  turned  off  so  that  the  total 
recording  capacity  of  the  oscilloscope  is  utilized  for  recording  the 
single  input  voltage.   This  signal  provides  the  information  necessary  for 
determining  the  pressure  history  of  the  test  element  during  the 
transient.   Experimental  runs  to  determine  the  test  element  temperature 
history  during  decompression  are  performed  by  recording  the  variations  in 
both  the  voltage  across  the  test  element  and  the  voltage  across  a 
standard  resistor.  Model  79  21  from  Ohmite,  with  the  oscilloscope 
triggered  by  the  latter  signal.   The  two  source  signals  are  fed  to  the 
voltage  bias  system,  where  the  voltages  Induced  by  the  application  of 
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power  at  initial  conditions  are  cancelled  by  applying  bucking  voltages. 
The  two  recorded  signals,  together  with  their  respective  bucking 
voltages,  provide  the  information  necessary  for  determining  the 
resistance  of  the  test  element  and  thus  its  temperature,  as  well  as  the 
power  delivered  to  the  test  element.   The  computer  code  used  in  the 
analysis  of  temperature  data  is  provided  in  Appendix  A,  and  the  computer 
code  used  in  the  analysis  of  pressure  data  is  provided  in  Appendix  B. 

4.5  Test  Element  Installation 

The  test  vessel  is  cleaned  using  acetone  and  then  rinsed  with  the 
test  fluid.   The  optical  ports  require  special  attention  during  cleaning 
to  keep  the  windows  clear  of  spots  and  also  free  from  scratches. 

The  test  element  is  soldered  to  the  brass  electrodes  of  the  test 
section.   The  entire  test  section  is  then  cleaned  with  acetone  to  remove 
excess  flux  or  other  contamination  arising  from  the  soldering  process. 
The  test  section  is  then  installed  in  the  pressure  vessel  and  the 
mounting  bolts  torqued  to  80  to  110  Nm.   A  small  voltage  is  applied 
across  the  test  element  in  such  a  way  that  the  wire  glows  cherry  red. 
The  air  heating  is  done  for  approximately  ten  minutes.   The  heating  of 
the  wire  tends  to  anneal  the  surface  of  the  wire  thus  leading  to 
consistent  surface  properties. 

After  the  Initial  preparation  of  the  test  element,  the  test  fluid  is 
poured  into  the  pressure  vessel.   The  wire  is  then  conditioned  by 
applying  a  small  voltage  across  it  until  boiling  initiated.   The  wire  is 
operated  at  this  voltage  for  approximately  ten  minutes,  then  the  voltage 
is  Increased  slightly.   This  is  repeated  until  the  voltage  across  the 
wire  reaches  a  maximum  of  10  V.   The  purpose  of  this  conditioning  is  to 
activate  the  cavities  on  the  wire  by  driving  out  any  trapped  gas. 
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Each  time  a  new  test  element  is  installed,  the  power  supply  must  be 
calibrated.   The  following  steps  are  used: 

1.  Set  V  and  V  equal  to  zero  by  grounding. 

2.  Adjust  R„  until  the  output  of  V  is  equal  to  zero, 

3.  Adjust  R„  until  the  output  of  V  is  equal  to  zero. 

4.  Adjust  R,  until  the  output  of  V   is  equal  to  zero, 

5.  Set  V  and  V,  equal  to  each  other  and  adjust  R,  so 

a      b   ^  -"A 

that  the  output  of  V  is  equal  to  zero. 

6.  Set  V,  to  zero  (leaving  V  alone),  then  adjust  R_  so 
that  the  output  of  V  is  equal  to  zero. 

7.  Set  V  equal  to  zero  and  re-adjust  R,  so  that  the 
output  of  V  is  again  equal  to  zero. 

8.  Set  V  equal  to  8.245  V  and  V  equal  to  4.588  V. 

9.  Adjust  R,  until  V  is  equal  to  10.00  V. 
10.   Adjust  R  until  V  is  equal  to  10.00  V. 

4,6  Procedures 

4.6,1  Preliminary  Procedures 

These  steps  are  followed  before  each  experimental  session: 

1.  The  test  fluid  is  heated  to  boiling  and  boiled  for  one  hour  to 
remove  dissolved  gases. 

2.  A  voltage  is  applied  to  the  test  element  until  vigorous  boiling 
has  commenced.   The  boiling  from  the  element  is  continued  for 
ten  minutes. 

3.  The  pressure  vessel  is  refilled  with  de-gassed,  near-boiling 
test  fluid,  and  all  valves  to  atmosphere  are  closed. 

4.  All  associated  electronic  equipment  is  allowed  to  warm  up  for 
the  hour  required  for  de-gassing. 
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5.   The  oscilloscope  Channel  B  trigger  level  is  adjusted  using  a 
Wavetek  Digital  VCG  Model  113. 

4.6.2  Experimental  Procedures 

Since  the  oscilloscope  has  only  two  input  channels  for  recording 
data,  pressure  data  and  temperature  data  must  be  recorded  in  separate 
experimental  runs.   However,  most  procedural  steps  are  the  same. 

These  steps  are  followed  to  set  up  each  experimental  run: 

1.  Vigorous  boiling  from  the  test  element  is  achieved  and 
continued  for  ten  minutes  at  atmospheric  pressure  and  the 
desired  bulk  fluid  subcooling.   This  allows  activation  of  the 
full  size  range  of  bubble  nucleation  sites. 

2.  The  test  element  is  allowed  to  cool  to  test  fluid  bulk 
temperature, 

3.  The  system  is  slowly  pressurized  to  the  desired  pressure  while 
holding  the  bulk  temperature  constant.   The  system  is 
maintained  at  this  condition  for  ten  minutes  to  allow  it  to 
stabilize. 

4.  Power  is  applied  to  the  test  element  until  the  desired  test 
element  temperature  is  achieved. 

5.  The  zero  settings  on  the  active  input  amplifiers  of  the 
oscilloscope  are  adjusted  and  the  trigger  time  zero  is  checked. 
These  settings  are  important  for  proper  triggering  of  the 
oscilloscope  and  for  use  in  the  computer  analysis  of  the  run. 

6.  Voltage  signals  from  the  standard  resistor  and  the  test  element 
are  fed  to  the  voltage  bias  system.   Bucking  voltages  are 
applied  such  that  the  signals  from  the  voltage  bias  system  to 
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the  oscilloscope  are  nominally  zero.   (This  step  is  performed 
only  if  recording  temperature  data.) 
If  the  boiling  initiation  time  is  not  desired  when  performing  a  run  to 
record  pressure  data,  steps  1,  2,  and  4  may  be  disregarded. 

After  the  run  has  been  set  up,  these  steps  are  followed  to  continue 
the  experiment: 

7.  Equipment  settings  are  checked  and  recorded  as  necessary,  and 
test  fluid  pressure  and  temperature  are  recorded.  The  gas 
supply  solenoid  valve  is  then  closed. 

8.  The  decompression  is  initiated  and  the  oscilloscope  is 
triggered  by  the  appropriate  input  signal. 

9.  The  recorded  data  are  stored  on  floppy  disc  for  later 
transmission  to  the  computer. 

10.  Manually  observed  boiling  initiation  time  is  recorded  (if 
measured) . 

11.  Steps  1  through  10  are  repeated  until  all  data  for  this 
combination  of  test  parameters  are  acquired. 

The  above  steps  are  repeated  for  other  test  parameter  combinations 
of  interest. 
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5.   RESULTS 

5.1  Boiling  Initiation 

The  initiation  of  boiling  at  a  surface  is  governed  by  the  size 
distribution  of  potential  nucleation  sites  and  the  fluid  temperature 
profile  adjacent  to  the  surface.   The  largest  site  will  be  the  first  to 
become  active  since  it  requires  the  least  superheat  for  boiling.   The 
mouth  diameter  of  the  largest  potentially  active  site  is  estimated  to  be 
0.35  urn.   This  value  was  obtained  from  photo-micrographs  of  the  surfaces 
of  similar  heating  elements  used  in  the  concurrent  study  of  pressure 
effects  on  boiling  initiation  during  transient  heating  (22) .   As 
discussed  in  Chapter  3,  the  activation  of  a  cavity  occurs  when  the 
superheat  has  reached  a  prescribed  magnitude  at  a  distance  from  the 
surface  comparable  to  the  radius  of  the  cavity  mouth. 

Natural  convection  determines  the  temperature  gradient  in  the  fluid 
adjacent  to  the  surface  before  decompression.   After  natural  convection 
is  fully  established,  the  thermal  boundary  layer  and  the  temperature 
profile  do  not  change  significantly  until  the  fluid  is  agitated,  either 
by  boiling  or  by  bulk  fluid  movement.   The  increase  in  superheat 
necessary  to  initiate  boiling  arises  from  the  decrease  in  the  bulk  fluid 
saturation  temperature  associated  with  the  decrease  in  system  pressure. 

The  temperature  of  the  test  element  should  behave  in  a  simple 
manner  when  a  constant  power  is  applied  and  natural  convection  from  the 
surface  to  the  fluid  is  fully  established.   The  temperature  would  be 
constant  before  a  pressure  transient  is  initiated,  and  would  remain  so 
during  the  transient  until  the  inception  of  boiling.   At  this  point  the 
test  element  temperature  would  begin  to  decrease.   Thus,  finding  this 
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point  of  time  in  the  temperature  history  and  the  pressure  at  this  time 
from  the  pressure  history  of  a  test  would  establish  the  conditions  for 
boiling  inception. 

In  this  study  of  boiling  initiation  during  pressure  transients, 
three  voltage  signals,  one  for  pressure  and  two  for  temperature,  were 
available  for  each  test.   However,  the  digital  oscilloscope  could  only 
record  two  of  these  simultaneously.   Since  both  temperature  signals  had 
to  be  recorded  for  quantitative  results,  the  nature  and  reproducibility 
of  the  decompression  event  had  to  be  investigated. 

The  system  was  prepared  for  experimental  runs  by  bulk  boiling  of 
the  water  in  the  pressure  vessel  for  one  hour  at  atmospheric  pressure, 
after  which  the  fluid  lost  through  evaporation  was  replaced  with 
degassed  water  at  near-boiling  temperature.   The  pressure  vessel  vent  to 
the  atmosphere  was  closed  after  bulk  boiling,  and  subsequent  system 
venting  to  the  atmosphere  was  accomplished  by  opening  the  gas  side  of 
the  pressurizer  to  atmospheric  pressure.   This  was  followed  by  vigorous 
boiling  from  the  platinum  test  element  for  ten  minutes.   Pressure  and 
temperature  runs  were  then  conducted  to  examine  the  effect  of  the 
pressure-temperature  history  on  boiling  initiation. 

5.2  Pressure  Runs 

The  program  used  to  analyze  the  pressure  runs  is  provided  in 
Appendix  B.   The  choice  of  amplifier  time  constant  and  the  effect  of 
system  temperature  on  the  pressure  signal  is  discussed  in  Appendix  C. 
Representative  results  from  the  program  are  provided  in  Appendix  D.   All 
pressure  runs  are  described  in  Table  5.1. 

The  reproducibility  of  the  decompression  was  of  primary  importance. 
In  all  runs  the  ambient  temperature  was  lOCC,  and  the  decompression 
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Table  5.1. 


Results  of  measurements  made  to  determine  the 
behavior  of  the  pressure  transients. 


Run 


(MPa) 


Drift  (V/s) 


Transfer  Function  (MPa/V) 


AP2-2 

0.377 

-0.005246 

AP2-6 

0.377 

-0.005246 

AP2-7 

0.515 

-0.0096 

AP2-8 

0.515 

-0.0096 

AP3-4 

0.515 

-0.0096 

AP3-5 

0.446 

-0.007801 

AP3-6 

0.446 

-0.007632 

AP4-2 

0.446 

-0.005812 

AP4-3 

0.584 

-0.006913 

AP4-7 

0.584 

-0.005800 

AP4-8 

0.584 

-0.005295 

AP7-5 

0.377 

-0.003981 

AP7-8 

0.377 

-0.003926 

AP8-6 

0.377 

-0.002703 

AP9-2 

0.377 

-0.004091 

AP9-8 

0.377 

-0.002093 

0.7942 
0.8654 

0.7798 
0.7908 
0.8131 

0.7631 
0.7705 
0.7748 


7631 
7771 


0.7690 

0.9952 
1.0297 
0.8536 
0.9432 
0.8248 
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ended  at  atmospheric  pressure.   The  pressure  histories  for 
decompressions  from  initial  pressures  of  0.377-0.584  MPa  are  illustrated 
in  Figs.  5.1  -  5.4,  where  it  is  seen  that  the  decompression  event  is 
reproducible  over  most  of  the  initial  pressure  drop. 

The  results  of  runs  AP2-7,  AP2-8,  and  AP3-4  require  special 
comment.   For  unknown  reasons,  analysis  of  the  raw  data  indicated 
negligible  drift  and  a  pressure  transducer  sensitivity  greater  than  the 
manufacturer's  specification.   Since  this  was  not  characteristic  of 
previous  experience,  a  drift  value  was  arbitrarily  specified  to  give  a 
transfer  function  reasonable  with  expectations.   However,  results  with 
this  data  should  be  viewed  with  caution. 

The  pressurizer  tended  to  back-fill  with  water  from  the  pressure 
vessel  in  the  course  of  the  repeated  fluid  temperature  cyclings 
experienced.   This  process  was  slow  enough  that  little  effect  on 
depressurization  was  observed  until  the  pressurizer  was  almost 
completely  full  of  water.   When  this  occurred  the  decompression  event 
became  erratic  and  unreliable.   When  approximately  two  liters  of  water 
were  removed  from  the  pressurizer,  allowing  a  larger  gas  volume  for 
pressurization,  the  rate  of  decompression  changed  significantly  (Fig. 
5.5).   However,  the  partial  removal  of  water  did  not  affect  the 
reproducibility  of  subsequent  pressure  runs,  as  shown  in  Figs.  5.6  and 
5.7. 

The  decompressions  displayed  near-exponential  behaviors  over  the 
reproducible  portions  of  the  transients.   Figures  5.8  and  5.9  illustrate 
this  behavior  for  representative  pressure  runs,  and  Table  5.2  lists 
their  approximate  reduction  periods.   These  periods  were  on  the  order  of 
4  s  before  the  pressurizer  was  partially  drained  and  on  the  order  of  6.6 
s  afterward. 


:.trf 


57 


D 

a. 


oc 


CO 
CO 


q:. 


.2  - 


TIME    (s) 


FIG.  5.1.   Experimental  pressure  results  for  decompression  from 

0.377  to  0.101  MPa  at  T  =  lOCC.   The  solid  line 
a 

represents  results  from  run  AP2-2.   The  dashed  line 
represents  results  from  run  AP2-6. 
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FIG.  5.2.   Experimental  pressure  results  for  decompression  from 

0.515  to  0.101  MPa  at  T  =  100°C.   The  dashed  line 
a 

represents  results  from  run  AP2-7,  the  solid  line 
represents  results  from  run  AP2-8,  and  the  dot-dashed 
line  represents  results  from  run  AP3-4. 


59 


a 

Q- 


UJ 
Ql 
ZD 
CO 

in 

ui 
cm 
a. 


FIG.  5.3.   Experimental  pressure  results  for  decompression  from 

0.446  to  0.101  MPa  at  T  =  100°C.   The  solid  line 
a 

represents  results  from  run  AP3-5,  the  dashed  line 
represents  results  from  run  AP3-6,  and  the  dot-dashed 
line  represents  results  from  run  AP4-2. 
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FIG.  5.4.   Experimental  pressure  results  for  decompression  from 

0.584  to  0.101  MPa  at  T  =  100°C.   The  solid  line 
a 

represents  results  from  run  AP4-3,  the  dashed  line 
represents  results  from  run  AP4-7,  and  the  dot-dashed 
line  represents  results  from  run  AP4-8. 
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FIG.  5.5.   Comparison  of  experimental  pressure  results  for 

decompression  from  0.377  to  0.101  MPa  at  T^  =  100°C 
before  and  after  partial  draining  of  the  pressurizer. 
The  solid  line  represents  run  AP2-2,  performed  before 
approximately  2  liters  of  water  were  removed  from  the 
pressurizer.   The  dashed  line  represents  run  AP8-8, 
performed  after  the  quantity  of  water  was  removed. 
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Table  5,2,   Approxmiate  reduction  periods  for  representative  pressure 
runs. 


Series 

Run 

tj  (s) 

t2  (s) 

p(tp  (MPa) 

p(t2)  (MPa) 

T  (s) 

1 

AP2-2 

0.56 

3.60 

0.3153 

0.1496 

4.078 

2 

AP2-8 

0.56 

3.60 

0.4300 

0,1921 

3.773 

3 

AP3-5 

0.56 

3.60 

0.3764 

0.1733 

3.920 

4 

APA-3 

0.56 

3.60 

0.4822 

0.2089 

3.634 

5 

AP7-5 

0.56 

5.52 

0.3449 

0.1736 

7.225 

6 

AP8-6 

0.48 

10.08 

0.3409 

0.1612 

6.622 

7 

AP9-2 

0.56 

5.52 

0.3396 

0.1587 

6.519 

8 

AP9-8 

0.56 

5.52 

0.3394 

0.1615 

6.680 
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5.3  Temperature  Runs 

The  preparation  of  the  system  for  a  series  of  runs  was  described 
earlier.   Additionally,  each  temperature  run  was  preceded  by  vigorous 
boiling  from  the  platinum  test  element  for  ten  minutes  at  atmospheric 
pressure  to  ensure  a  consistent  range  of  cavity  sizes  before 
pressurlzatlon. 

In  the  series  of  tests  presented,  the  system  was  pressurized  to  a 

maximum  ambient  pressure,  p     ,  while  the  bulk  water  temperature  was 

held  steady  at  100°C.   Power  was  then  applied  to  the  test  element  to 

elevate  its  temperature  to  160°C.   The  decompression  event  was  then 

Initiated  and  a  manual  measurement  of  boiling  initiation  time  t.  was 

b 

obtained.  The  program  used  to  analyze  the  temperature  runs  is  provided 
in  Appendix  A.  Representative  results  from  the  program  are  provided  in 
Appendix  E.   All  series  are  described  in  Table  5.3. 

In  series  1  the  system  was  pressurized  to  a  maximum  pressure  of 
0.377  MPa.   Figure  5.10  illustrates  series  1,  where  It  is  observed  that 
the  temperature  of  the  test  element  began  to  decrease  almost 
immediately,  while  the  measured  boiling  initiation  time  occurred  after 
the  temperature  had  dropped  an  appreciable  amount. 

In  series  2  the  system  was  pressurized  to  0.515  MPa.   The  results 
of  this  series  are  shown  in  Fig.  5.11.   The  temperature  of  the  test 
element  remained  constant  for  approximately  two  seconds,  and  the 
measured  boiling  times  were  close  to  the  times  at  which  the  temperatures 
begin  to  decrease.   Also  apparent  is  the  variability  of  the  heater 
temperature  recovery  to  a  steady-state  value. 

Figures  5.12  and  5.13  illustrate  series  3  and  series  4,  in  which 
the  system  was  pressurized  to  0.446  and  0.584  MPa,  respectively.   Again, 
the  heater  temperature  remained  constant  for  a  period  of  time  after 
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Table  5.3.   Experimental  conditions  for  measurements  to  determine  the 

effects  of  pressure-temperature  history  on  boiling  initiation. 
Initial  conditions:   wire  temperature  T  =  150°C,  ambient 
temperature  T  =  100°C.  " 


Series 


Run 


Duration  (s) 


p      (MPa) 
a, max 


p  (t=0)  (MPa)    q  (MW/m^) 


AP2-3 

20 

0.377 

0.377 

0.427 

AP2-4 

20 

0.377 

0.377 

0.430 

AP2-5 

20 

0.377 

0.377 

0.432 

AP3-1 

20 

0.515 

0.515 

0.424 

AP3-2 

20 

0.515 

0.515 

0.424 

AP3-3 

20 

0.515 

0.515 

0.437 

AP3-7 

20 

0.446 

0.446 

0.419 

AP3-8 

20 

0.446 

0.446 

0.419 

AP4-1 

20 

0.446 

0.446 

0.419 

AP4-4 

20 

0.584 

0.584 

0.408 

AP4-5 

20 

0.584 

0.584 

0.416 

AP4-6 

20 

0.584 

0.584 

0.416 

AP7-2 

20 

0.377 

0.377 

0.407 

AP7-3 

20 

0.377 

0.377 

0.419 

AP7-4 

20 

0.377 

0.377 

0.419 

AP7-6 

20 

0.377 

0.377 

0.419 

AP7-7 

20 

0.377 

0.377 

0.419 

AP8-4 

20 

1.480 

0.377 

0.419 

AP8-5 

20 

>1.480 

0.377 

0.414 

AP8-7 

20 

1.480 

0.377 

0.418 

AP8-8 

100 

1.480 

0.377 

0.424 

AP9-4 

200 

0.791 

0.377 

0.418 

AP9-6 

200 

0.791 

0.377 

0.418 
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initiation  of  the  pressure  transient,  and  measured  boiling  times  were 

close  to  the  points  where  the  temperatures  began  to  decrease.   Figure 

5.13  also  shows  particularly  well  the  variable  nature  of  the  heater 

temperature  recovery  after  boiling  had  begun. 

Series  5-8  were  performed  after  removal  of  a  quantity  of  water  from 

the  pressurizer.   The  rates  of  decompression  in  these  runs  were  less 

than  those  for  series  1-4  as  discussed  earlier.   Additionally,  the 

procedure  was  varied  such  that  the  system  was  subjected  to  a  maximum 

pressure  p     ,  then  slowly  returned  to  a  starting  pressure  p  (t=0), 
a,  max  3. 

the  pressure  at  the  initiation  of  the  pressure  transient. 

In  series  5  the  maximum  pressure  was  0.377  MPa.   Measured  boiling 
initiation  times  for  all  runs  in  this  series  except  for  run  AP7-6  were 
well  after  the  heater  temperature  had  begun  to  decrease.   As  in  series 
1,  the  heater  temperature  began  to  decrease  almost  immediately  after 
initiation  of  the  pressure  transient.   Figure  5.14  illustrates  this 
series. 

In  series  6  and  series  7  the  maximum  pressure  applied  was  1.480 
MPa.   The  exception  was  run  AP8-5,  which  experienced  an  unknown  maximum 
pressure  greater  than  1.480  MPa.   In  these  runs  the  measured  boiling 
time  corresponds  to  the  time  of  the  first  audible  bubble  collapse. 
Bubbles  did  not  appear  on  the  visible  part  of  the  test  element  until 
after  the  system  had  returned  to  atmospheric  pressure.   In  series  6 
(Fig.  5.15),  the  heater  temperature  began  to  decrease  shortly  after 
initiation  of  the  pressure  transient,  but  the  decrease  was  very  gradual 
until  the  end  of  the  pressure  transient.   Heater  temperature 
fluctuations  indicative  of  substantial  nucleate  boiling  then  appeared, 
but  the  temperature  decreased  minimally  compared  to  runs  with  less 
extreme  overpressures . 
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This  behavior  was  also  observed  In  series  7,  in  which  runs  of  20  s  and 
100  s  duration  were  performed.   The  temperature  behavior  during  the 
first  20  s  of  the  event  is  shown  in  Fig.  5.16(a).   Even  100  s  after  the 
initiation  of  the  pressure  transient  [Run  AP8-8  and  Fig.  5.16(b)]  the 
heater  temperature  had  only  decreased  about  7  K,  as  opposed  to  10-12  K 
drops  in  20  s  seen  in  series  5,  for  example. 

The  long-term  behavior  of  the  heater  temperature  was  of  specific 
interest  in  series  8.   The  maximum  pressure  in  this  series  was  0.791 
MPa,  and  the  duration  of  the  runs  was  200  s.   Again,  no  bubbles  were 
seen  until  after  the  system  had  returned  to  atmospheric  pressure. 
Measured  boiling  initiation  times  were  before  substantial  nucleate 
boiling  was  indicated  by  the  temperature  analysis.   The  minimal  decrease 
in  the  heater  temperature  during  the  first  20  s  of  the  event  is 
illustrated  in  Fig.  5.17(a).   After  the  initial  drop,  the  temperature 
reached  a  plateau  before  experiencing  a  more  sustained  decrease.   The 
heater  temperature  approached  the  steady-state  boiling  temperature  at 
the  same  superficial  heat  flux  by  approximately  200  s.   Series  8  and  the 
corresponding  steady-state  boiling  temperature  traces  are  shown  in  Fig. 
5.17(b). 

If  the  point  at  which  the  heater  temperature  starts  to  decrease  is 
taken  as  the  incipient  boiling  point,  conditions  at  boiling  initiation 
can  be  determined  from  pressure  and  temperature  values  at  this  time. 
Normally,  this  estimated  time  should  be  equivalent  to  the  measured 
boiling  Initiation  time,  with  any  difference  accounted  for  by  resolution 
of  the  data  analysis,  human  error,  and  reaction  time  in  operating  the 
chronograph.   However,  a  significant  discrepancy  could  exist  between 
these  two  times.   If  the  resolution  limit  of  the  analyzed  data  Is 
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approached,  the  measured  value  should  be  selected  as  the  true  boiling 
initiation  time.   If  boiling  initiates  on  a  part  of  the  test  element 
outside  the  visual  range  and  the  ambient  noise  level  precludes  audio 
detection,  the  measured  time  would  be  late  and  the  estimated  value 
should  be  selected  as  the  true  boiling  initiation  time.   Additionally, 
although  the  partial  pressure  of  noncondensible  gas  is  assumed  to  be 
negligible,  an  isolated  cavity  could  contain  a  bubble  of  noncondensible 
gas.   Premature  non-boiling  bubble  formation  could  occur  at  this  site. 
The  frequency  of  bubble  formation  would  not  display  the  periodicity  of 
true  nucleate  boiling.   Although  conditions  for  boiling  initiation  would 
not  have  been  achieved  according  to  the  criteria  of  this  work,  the 
chronograph  would  have  been  stopped  at  the  appearance  of  the  first 
bubble.   Again,  the  estimated  boiling  initiation  time  should  be 
selected.   Following  this  prescription  for  determining  true  boiling 
initiation  times,  the  measured  and  true  boiling  initiation  times  and 
conditions  at  boiling  initiation  are  given  in  Table  5.4  for  all  series. 
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6.   DISCUSSION  OF  RESULTS 

6.1  Pressure  Measurement 

After  the  initiation  of  decompression  the  pressure  vessel 
experienced  a  near-exponential  reduction  in  pressure.   Pressure 
transients  were  generally  reproducible  to  times  well  beyond  measured 
boiling  initiation  times.   The  spikes  that  occasionally  appeared  in  the 
measurements  are  believed  to  be  due  to  spurious  electronic  noise.   The 
pressure  measurement  system  sensitivity  was  less  than  expected  due  to 
installation  constraints  (see  Appendix  C) .   However,  the  system 
performed  consistently  except  in  the  decompressions  from  0.515  MPa 
(series  2) .   The  reasons  for  the  unusual  behavior  in  those  runs  are 
unknown;  therefore  the  results  of  series  2  are  not  considered  in  the 
following  discussion, 

6.2  Estimation  of  Boiling  Initiation  Times 

Consider  a  heater  experiencing  a  constant  power  application.   In 
the  case  where  heat  transfer  from  the  heater  to  a  surrounding  fluid  is 
by  natural  convection,  the  inception  of  nucleate  boiling  will  enhance 
the  heat  transfer  and  cause  the  heater  surface  temperature  to  decrease. 
Estimated  boiling  initiation  times  were  obtained  by  determining  the 
times  at  which  the  test  heater  temperature  began  to  show  a  substantial 
decrease.   These  times  provided  an  independent  means  of  confirming 
measured  boiling  initiation  times.   The  accuracy  of  the  estimated  times 
was  dependent  on  the  techniques  used  to  obtain  and  analyze  the  raw 
voltage  data.   Since  the  variations  in  the  voltages  were  expected  to  be 
small  compared  to  the  initial  applied  voltages,  the  measurement  system 
was  set  up  to  measure  only  the  variations.   The  analysis  of  the  data 
required  averaging  the  data  over  several  points  to  minimize  masking  of 
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changes  in  true  voltage  by  electrical  noise.   However,  the  data 
reduction  resulted  in  coarse  time  steps  and  thus  crude  estimations  of 
the  time  at  which  the  heater  temperature  began  to  decrease.   Although 
the  two  methods  of  determining  boiling  initiation  times  should  have 
produced  the  same  values,  this  was  observed  only  in  series  3  and  series 
4,  where  the  inception  of  boiling  was  accompanied  by  a  marked  change  in 
heater  temperature.   In  these  two  series,  the  true  boiling  initiation 
times  were  the  measured  times  adjusted  according  to  the  temperature 
histories.   These  adjustments  were  reasonable  considering  reaction  times 
and  human  error.   In  the  other  series,  however,  experimental  runs 
generally  lacked  a  sharp  change  in  heater  temperature  at  boiling 
initiation,  and  the  combination  of  test  conditions  helped  to  produce 
major  differences  in  estimated  and  measured  times.   In  series  1  the 
temperature  began  to  decrease  almost  immediately,  so  true  boiling 
initiation  times  were  the  first  times  for  which  temperature  data  were 
calculated.   In  series  5-8  poor  data  resolution  combined  with  relatively 
slow  pressure  transients  required  specification  of  a  range  of  boiling 
initiation  times.   Table  6.1  summarizes  the  conditions  at  boiling 
initiation  for  the  series  considered. 

6.3  Effects  of  Pressure-Temperature  History  on  Boiling  Initiation 

Runs  in  each  series  had  been  preceded  by  vigorous  boiling  from  the 
test  element  for  10  minutes  at  atmospheric  pressure  and  bulk  fluid 
temperature  T  =  100°C.   The  pressure  was  raised  to  the  maximum  value 
for  each  run  before  power  was  applied  to  the  test  heater  to  elevate  its 
temperature. 

In  series  1,  3,  and  4  the  decompression  began  at  the  maximum 
pressure,  and  reduction  periods  were  on  the  order  of  4  s.   The  boiling 
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Table  6.1   Summary  of  results  of  measurements  made  to  determine  the  effects  of 

pressure-temperature  history  on  boiling  initiation.   In  all  cases,  the 
heater  is  a  platinum  wire  of  0.25  mm  diameter  and  9.6  cm  length. 
Degassing  was  carried  out  at  atmospheric  pressure.   Initial  conditions: 
T  =  160°C,  T  =  100°C. 


Series 

p               (MPa) 
a, max 

p^(t=0)    (MPa) 

%^ 

s) 

Pb    ' 

[MPa) 

b      s 

(K) 

1 

0.377 

0.377 

0.08  - 

0.40 

0.365 

-   0.330 

19.7   - 

23.2 

3 

0.446 

0.446 

1.84  - 

2.00 

0.269 

-   0.259 

30.2  - 

31.4 

4 

0.584 

0.584 

3.12  - 

3.92 

0.236 

-   0.194 

34.4  - 

40.7 

5 

0.377 

0.377 

1.0     - 

2.5 

0.321 

-   0.260 

24.1   - 

31.3 

6 

1.48 

0.377 

8.88  - 

12.24 

0.110 

-   0.101 

57.7   - 

60.0 

7 

1.48 

0.377 

7.6      - 

9.2 

0.127 

-  0.114 

53.6   - 

56.7 

8 

0.791 

0.377 

6.4     - 

8.0 

0.144 

-  0.121 

49.9  - 

55.0 

S6 

initiation  times  increased  with  Increasing  initial  pressure,  and  no 
overlapping  of  times  was  observed  for  the  different  series.   Increasing 
initial  pressure  also  resulted  in  decreasing  boiling  initiation 
pressures  and  increasing  superheats  at  boiling  inception. 

In  series  1,  Fable's  hypothesis  gives  the  largest  available 
nucleation  sites  radii  of  r  =  0.132  um.   This  would  lead  to  a  boiling 
initiation  pressure  of  0.330  MPa  and  a  boiling  initiation  superheat  of 
23.2  K.   In  series  3,  r^  =  0.106  ym,  leading  to  a  boiling  initiation 
pressure  of  0.258  MPa  and  a  boiling  initiation  superheat  of  31.5  K. 
Both  of  these  series  are  in  excellent  agreement  with  the  Fable  model  if 
the  estimated  boiling  inception  times  are  used.   However,  the  model 
predicts  a  boiling  initiation  pressure  of  0.114  MPa  for  the  cavity  size 
r^  =  0.075  um  of  series  4.   This  pressure  is  substantially  below  that 
observed  in  this  series. 

In  series  5-8  the  pressure  transient  was  initiated  from  0.377  MPa 
after  the  initial  maximum  pressure  was  applied.   Reduction  periods  were 
on  the  order  of  6.6  s.   Boiling  initiation  times  increased  with 
increasing  initial  pressure,  but  the  range  of  times  tended  to  overlap 
for  maximum  pressures  of  0.791  and  1.48  MPa. 

Series  5  repeated  the  conditions  of  series  1  at  a  slower 
decompression  rate.   Although  superheats  were  higher  and  pressures  were 
lower  at  the  inception  of  boiling  in  series  5,  the  results  also  support 
Fable's  model.   Series  6-8  experienced  such  extremes  in  pressure  that 
Fable's  model  predicted  no  boiling  for  the  given  heater  temperature. 
The  model  predicted  a  maximum  cavity  radius  of  0.026  um  for  series  6  and 
series  7,  which  would  require  a  heater  temperature  of  192. 3°C  to 
initiate  boiling.   Series  8  would  require  a  heater  temperature  of 
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169. T'C  to  initiate  boiling  from  a  predicted  maximum  cavity  radius  of 
0.053  pm.   However,  the  heater  temperature  history  indicated  that 
boiling  began  near  the  end  of  decompression,  and  after  an  additional 
period  of  time  sites  were  seen  to  commence  nucleating  on  the  visible 
part  of  the  test  element. 

6.4  Approach  to  Steady-State  Boiling 

Even  though  Fable's  model  greatly  underpredicted  boiling  initiation 
pressures  (and  overpredicted  boiling  initiation  superheats)  in  series  4 
and  series  6-8,  the  pressure-temperature  history  did  affect  the 
reactivation  of  nucleation  sites  as  the  runs  progressed. 

The  models  used  to  predict  boiling  conditions  based  on  cavity  size 
can  also  explain  the  behavior  of  the  heater  temperature  as  steady-state 
nucleate  boiling  is  approached.   Consider  Eq .  (2.7), 

Pv^V  -  Pv"s)  =  2a(y/r^.  (2.7) 

Although  surface  tension  decreases  with  increasing  temperature,  vapor 
pressure  is  more  sensitive  to  temperature  variations  and  hence  dominates 
the  temperature  dependence  of  Eq.  (2.7).   It  is  apparent  that  test 
element  temperature  results  would  be  very  sensitive  to  whether  or  not 
one  or  two  cavities  activate  consistently  in  separate  runs. 

Since  T  remains  constant  until  the  Inception  of  boiling,  superheat 
increases  with  decreasing  pressure  until  conditions  for  boiling  are 
achieved.   Once  boiling  begins,  the  heater  cools,  with  the  amount  of 
cooling  determined  by  the  size  and  density  of  active  sites.   During  the 
pressure  transient,  the  cooling  tends  to  offset  the  decrease  in  p  (T  ) , 
thus  restricting  the  activation  of  other  nucleation  sites. 
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After  the  decompression  ends,  a  larger  site  rendered  inactive  by 
the  pressure-temperature  history  experienced  prior  to  the  pressure 
transient  could  be  activated  in  two  ways.   Vapor  microbubbles  created  in 
the  collapse  of  macroscopic  bubbles  in  subcooled  liquid  could  migrate 
from  an  adjacent  active  site  and  seed  the  cavity  (33-36).   If  the  cavity 
has  vapor  trapped  at  its  base,  activation  could  be  accomplished  either 
by  vapor  migration  or  by  achievement  of  a  sufficiently  large  local 
superheat  to  overcome  the  effects  of  pre-pressurization.   Once  the 
barrier  to  nucleation  were  breached,  the  cavity  would  remain  active. 

The  temperature  of  the  platinum  wire  test  element  was  calculated 
using  its  resistance  at  0°C  (see  Section  4.2).   Based  on  measurements  of 
the  resistance  at  100°C,  the  wire  resistance  at  0°C  was  R  =  0.1961  ± 
0.0021.   This  led  to  an  uncertainty  of  ±2.9  K  in  the  value  of  the 
calculated  wire  temperature  T  .   This  was  a  systematic  error  that  adds 
to  the  uncertainty  in  boiling  initiation  conditions  beyond  that  reported 
in  Table  6.1. 

The  system  pressure  during  the  decompression  was  calculated  using 
the  recorded  value  of  the  pressure  at  the  initiation  of  the  event.   The 
error  in  this  measurement,  including  the  ±0.25%  accuracy  of  the  pressure 
gauge,  was  estimated  to  be  ±3.4  kPa.   In  addition,  the  pressure 
transducer  had  a  linearity  of  ±1.0%  (based  on  the  manufacturer's 
performance  data) ,  and  the  dual  mode  amplifier  had  a  range  accuracy  of 
±1%.   The  combination  of  these  sources  of  error  led  to  an  uncertainty  in 
the  pressure  measurement  of  approximately  ±8-9  kPa,  dependent  upon  the 
magnitude  of  the  pressure  drop  recorded.   This  error  also  adds  to  the 
uncertainty  in  boiling  conditions  beyond  that  reported  in  Table  6.1. 
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6.5  Conclusions 

Although  this  study  did  not  allow  a  quantitative  comparison  of 
results  to  prediction  based  on  Fable's  theory,  the  results  did  show  a 
clear  effect  of  the  pressure-temperature  history  on  the  initiation  of 
boiling  during  pressure  transients.   Increasing  the  initial  pressure 
tended  to  delay  the  inception  of  boiling,  showing  that  lower  pressures 
and  higher  superheats  were  necessary  to  initiate  boiling.   Since  smaller 
cavities  require  greater  superheat  to  nucleate,  this  is  consistent  with 
predictions  that  the  size  of  the  largest  potential  nucleation  site 
becomes  progressively  smaller  as  greater  overpressures  are  applied. 
That  boiling  occurred  in  cases  where  none  was  predicted  indicates  that  a 
simple  time- independent  theory  using  constant  contact  angles  is 
inadequate. 
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7.   SUGGESTIONS  FOR  FURTHER  STUDY 

A  more  comprehensive  study  of  the  pressure  and  temperature  effects 
on  boiling  initiation  during  pressure  transients  needs  to  be  performed, 
with  a  wide  range  of  conditions  and  different  heater  geometries. 
Certain  changes  are  recommended  to  improve  the  validity  and  efficiency 
of  such  a  study. 

The  test  heater  power  supply  should  be  upgraded  with  additional 
deep-cycle  storage  batteries.   This  would  improve  stability  and  allow 
application  of  greater  power  levels.   More  stable  power  supplies  to 
provide  biasing  voltages  are  also  necessary  to  enhance  the  resolution  of 
the  small  voltage  variations  associated  with  this  study. 

An  electro-acoustical  device  should  be  used  in  an  independent 
method  for  detecting  the  onset  of  nucleate  boiling.   The  visual 
technique  employed  in  this  study  was  hampered  by  a  limited  field  of 
vision  and  reaction  times.   Visual  detection  also  becomes  increasingly 
unreliable  as  pressure  reduction  periods  decrease. 

A  much  more  sensitive  pressure  transducer  should  be  installed. 
This  would  greatly  improve  the  response  of  the  pressure  measurement 
system  and  minimize  amplifier  drift  problems. 

A  four-channel  recording  digital  oscilloscope  and  a  captive 
computer  system  are  also  recommended.   With  a  four-channel  oscilloscope, 
all  pertinent  data  for  each  experimental  run  could  be  recorded,  and  the 
reproducibility  of  the  pressure  transient  would  not  be  a  major  concern. 
A  captive  computer  system  would  greatly  Improve  the  efficiency  of  the 
study,  allowing  each  test  to  be  analyzed  during  preparation  for  a 
subsequent  experiment.   Additionally,  near  immediate  feedback  of  the 
effects  of  changes  in  test  conditions  would  be  available. 
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Finally,  a  study  of  decompression  with  heater  temperature  held 
constant  needs  to  be  performed,  as  well  as  an  Investigation  of  the 
effect  of  temperature  on  advancing  and  receding  contact  angles.   These 
two  subjects  need  to  be  addressed  in  order  to  develop  a  general  model 
for  boiling  behavior  during  transient  condition. 
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APPENDIX  A 
Program  BOIL:   Temperature  Data  Analysis 

This  program  was  written  in  Hewlett-Packard  BASIC,  Version  2.0,  and 
run  on  the  HP  9816  microcomputer.   The  program  was  originally  written 
for  analysis  of  power  transient  data  from  a  concurrent  study  (22),  and 
this  capability  is  retained  in  the  current  version.   In  fact,  the 
program  parameters  default  to  the  power  transient  values,  and  changes  In 
these  values  arising  from  choosing  the  pressure  transient  option  are 
contained  in  the  voltage  input  data  set. 

The  voltage  data  for  an  experimental  run  are  measured  by  a  Nlcolet 
Explorer  III  digital  oscilloscope.   The  voltage  input  is  converted  to  a 
corresponding  channel  number  and  stored.   The  channel  number  is  obtained 
by  dividing  the  sampled  voltage  by  the  conversion  factor 

\orm  '   '^°1«38'=  Range}/2000     Volts/channel. 

There  are  4096  voltage  channels  (-2048  to  +2047),  which  allows  a  full 
scale  voltage  interval  greater  than  -VR  to  +VR,  where  VR  is  the 
oscilloscope  voltage  range  setting.   There  are  also  4096  memory 
addresses  available.   Channel  numbers  are  stored  sequentially  at  a  time 
interval  determined  by  the  Time  Per  Point  switch.   If  both  input 
channels  A  and  B  are  active,  the  signals  are  stored  in  the  sequence 
ABABAB,  and  the  time  interval  defines  the  time  between  each  sampling  of 
the  pair  AB.   Thus,  the  channel  number  defines  the  voltage  and  the 
memory  address  defines  the  time  at  which  the  voltage  was  sampled. 

Voltage  data  is  originally  input  to  the  program  from  the 
oscilloscope.   The  data  Is  also  saved  on  a  disk  for  reference.   The 
information  required  from  each  run  depends  on  the  case  considered,  i.e.. 
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transient  power  or  transient  pressure.   The  former  is  always  performed 
in  Normal  Trigger  mode  on  the  oscilloscope,  so  the  first  datum  is  at 
time  zero,  while  the  latter  is  always  performed  in  the  Cursor  Trigger 
mode.   In  this  configuration  the  first  datum  time  and  the  event 
initiation  time  must  be  recorded,  as  well  as  certain  equipment  settings. 
These  values  are  stored  as  the  first  nine  points  in  the  data  file  when 
it  is  saved  on  the  disk. 

The  program  separates  the  raw  voltage  data  into  the  two  input 
signals,  converts  the  data  to  true  voltages,  and  then  averages  the  data. 
The  voltage  is  found  from  the  expression 


V 

V  =  (Channel  number)  "°"°  -  Vcz  +  V 
Amp  o 


where  Amp  is  the  signal  amplification,  Vcz  is  the  zero  input  voltage 
level  for  the  oscilloscope,  and  V  is  the  biasing  voltage.   The  biasing 
voltage  is  subtracted  from  the  signal  before  it  is  amplified  and  sent  to 
the  oscilloscope  to  permit  greater  precision  in  the  measurement  of  the 
small  voltage  changes  associated  with  this  case.   The  voltages  are 
smoothed  over  n  points,  and  the  averaged  voltage  is  assigned  the 
averaged  time  value.   Smoothing  intervals  are  exclusive,  i.e.,  no  datum 
is  used  more  than  once  in  the  averaging  process. 

Program  output  is  available  in  an  ASCII  file  for  hardcopy  printout 
or  in  binary  data  files  for  plotting  the  results. 
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and  b)  CURSOR  TRIGGER  OF  0-SCOPE. 
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KSU   SEP  1984 


OPTION  BASE  1 

CON  /Array/  Raii'4096) ,A(512) ,B!512) ,EI512) ,2 1512, II 

CON  /Const/  Alph,fleta,Gbeg,l(cnd 

CON  /Datui/  Di,Len,Ra,R5,Ta,Tau,T0,Va0,Vb0,Vc:l,«C!2,Vnori,PJ[35] 

CON  /Stat/  Aip,node,Nav.NblDcl(,Si9«a 

Dllt  Title«[80i 

PRINT  CHR«(i:i 

INPUT  'RESET  GRAPHICS  (V/NI  ?', Reset* 

INPUT  'HHAI  IS  THE  RUN  NUMBER,  e.g.,  RUN  ?."?", Title* 

Title«="RUN  NUKBER  •ITitlel 

PRINT  Title! 

NAIN  PROGRAH 

CALL  Para« 

CALL  InptlN«) 
Tscale=Tsul2000 

IF  Reset»="V  THEN  CALL  6raf  setiTscale) 
PRINT  Titlel,'  ','  ', "SOURCE  DATA  FILE;  ';NJ 
PRINT 
PRINT  'AVERAGING  DATA  >  "iNavi"  POINT  AVERAGING' 

CALL  Avg 
PRINT  "CALCULATING  INITIAL  HIRE  RESISTANCE' 

CALL  Resis 
PRINT  'CALCULATING  PROPERTIES" 

CALL  Props 

CALL  OutpiN*. Title*) 
END 


SUB  Par  an 
CON  /Const/  Alph.BetajSbeo.Kcnd 

CON  /DatuD/  Di.L»n,Ra,R5,Ia,Tau,TO,VaO.VbO,Vczl,Vc:2,Vnor«,F« 
CON  /Stat/  A«p,Node,Nav,Nblocl(,Siq«a 

'   DAIA  FOR  THE  HIRE,  THE  FLUID,  AND  THE  SETTINGS  ON  THE  SCOPE 
I   DEFAULT  FOR  CONSTANT  PRESSURE  RUNS 


R!=. 2394*. 00203 

Dl=2.5E-4 

Len=?.6E-2 

Alph=. 00392 

Beta=-5.5E-7 

lcnd=.6740 


STANDARD  RESISTANCE 
HIRE  DIAMETER 
NIRE  LENGTH 


BEAN  VALUE  70  TO  200  C  IS  ,674  HATER  COND) 


98 


590  Vnori-.002                             !  4V  RSN6E  ON  O-SCOPE 

600  Vi:zl=-4.                                 '  ZERO  INPUT  VOLTAGE  LEVEL 

610  Vc:2M. 

620  VjO=0.               !  BUCKING  VOLTAGE 

630  VbO=0. 

640  I 

650  Sbe5=0. 

660  l1ode=l 

670  fl»p=l. 

680  Nav=16 

690  SUBEND 

700 

710 

720 

730  SUB  Inpt(N«) 

740  OPTION  BASE  I 

750  con  /Array/  R3iii4096),A(5I2),Bi5I2),Ei512l,ZI512,4) 

760  con  /Const/  Alph,Beta,Bbeo,Xcnd 

770  con  /Batui/  l)i,Len,Ra,P5,Ia,Iau.IO,VaO,VbO,Vczl,Vc:2,Vnori,P$ 

780  con  /Stat/  Ann, node, Nav.NblockiSigiia 

790  Din  Ra»«(4096f[51 

800  I 

810    INPUT  "DATA  SOURCE SCOPE  =  1,  BOAT  FILE  =  2", It 

820  IF  It=l  THEN 

830  PRINT 

840  PRINT  'TRIGGER  nODE NORMAL  =  1,  CURSOR  =  T 

850  INPUT  'INORnAL  if  constant  pressure  run)      ".node 

360  END  IF 

870  CALL  Reai)er(Tt,N$,Ra»l(lll 

880  IF  Ra««(ll='CURSR"  THEN  node=2 

890  IF  node(>2  THEN  Raii(l)=VAL(Ra»«iI)l 

900  FOR  1=2  TO  4096 

910  Raul  I)=VAL!Rai(«(  III 

920  NEH  I 

930  IF  node=l  THEN  1040 

940  Vnori=Raii(2l 

950  ChO=Ra«(3l 

960  Taii=Ra«l4) 

970  Vczl=RaK(51/1000 

980  Vc!2=Rai((6)/1000 

990  A«p=Raii(7l 

1000  VaO=Rait(8) 

1010  Vb0=-Rai<(9l 

1020  Reserv=16 

1030  P«="" 

1040  BEEP  2197. 26.. 3 

1050  INPUT  'ENTER  ANBIENT  TEnPERAIURE  ICI'.Ta 

1060  INPUT  "ENTER  INITIAL  NIRE  lENPERATURE  (C)'.TO 

1070  IF  node=2  THEN  GOTO  1110 

1080  INPUT  'ENTER  AHBIENT  PRESSURE,  e.g.  101.4.475.8,1397.8  (kPa)',P« 

1090  INPUT  'ENTER  INITIAL  HEAT  FLUI  ON  THE  HIRE  IN/i-2)',Bbeg 

1100  INPUT  "ENTER  SAnPLING  TINE  INCREHENT,  £.6.  lE-2  FOR  20  SEC  RUN'.Tau 

1110  INPUT  "ENTER  [4,8,16!  POINT  RUNNING  AVERAGE  (DEFAULT  =  lA)',Nav 

1120  IF  «ode=l  THEN  1250 

1130  INPUT  'ENTER  (lAIIHUn  PRESSURE', Pnas 

1140  INPUT  'ENTER  STARTING  PRESSURE  FOR  DECOnPRESSION',Pstart 

1150  INPUT  'UNITS  ON  VALUES  JUST  ENTERED  K  =  kPa,  P  =  psio',Unitt 

1160  IF  Unit»='r  THEN  1210 

1170  Corv=(tP»ax/14.696H-l)tI01.325 

1180  Piax'DR0UND(Conv,5l 

1190  Conv=((P5tart/14.696)*lin01.325 

1200  Pstart=DR0UND(Cony,5) 

1210  IF  Pstart<>P«a«  THEN  Pt="na»  "SVALIIPnasli",  " 

1220  P«=PM'DecDiipre55  (ro«  "MALJIPstartl 

1230  Neg=INIl(ChO-l-Re5erv)/Nav/2)        !  ALLONS  CALCULATIONS  FOR  TinE<0 

1240  Ra»(10)=ChC-NavlNeql2                     I  START  CHANNEL 

1250  Nch=(4094-iRai<(10)-l)l/2 

1260  NbIocl!=lNT(Nch/Nav) 

1270  PRINT  CHR«(12) 

1230  PRINT  Rait«ll),Raii(2),Rait(3),Rai.(4),Rai((5i,Rai((6),Raiil7),Raii(8) 

1290  PRINT 

1300  SUBEND 
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1310  ! 

1320  ! 

133(1  ' 

1340  SUB  Reader  iTt.Nane*,  Raul  ill) 

1350  OPTION  BASE  1 

1360  con  /Stat/  Hip, node, Nav,Nblock,Sigiia 

137(1  ! 

1330  PRINT 

1370  PRINT  'PLEflSE  INSERT  DATA  DISC  INTO  DISC  DRIVE  0' 

1400  INPUT  'PLEASE  ENTER  A  UNIQUE  NAME  FOR  THE  BOAT  fILE.",Na«e« 

1414  DISP  -NORKING,  PLEASE  NAIT." 

1420  NASS  STORAGE  IS  ":HP8:90I, 700,0" 

H30  IF  Tt=2  THEN  1650 

1440  ! 

1450  ASSIGN  »Scope  TO  S 

1440  CONTROL  9,3;7600 

1470  ASSIGN  eScopeiPORNAT  ON 

1480  CONTROL  9, 5; 3 

1490  CONTROL  9,4|2 

1500  OUTPUT  SScopejCHRKl); 

1510  OUTPUT  ?ScDpeiCHR«(69l!CHR«(48)iCHR«(68l!CHR«(49l!CHRI(6fl)!CHR$(48): 

1520  OUTPUT  »Sc!!pe;CHR$(79)iCHR«(52l;CHR»(48);CHR«(57)|CHR«(54) 

1530  OUTPUT  »ScopeiCHR«(2l: 

1540  ENTER  »Scope  USING  ■5A,l(,X";Raii«(ll 

1550  CONTROL  9,5:0 

ISiO  ASSIGN  eScope  TO  t 

1570  IF  Mode=2  THEN  CALL  Scope  set (Raw* 1 1) I 

1580  DIBP  'WORKING,  PLEASE  NAIT.' 

1590  CREATE  BOAT  Naie«,4096,5 

two  ASSIGN  $Path  TO  Naiet 

Ul«  OUTPUT  »Path  USING  'SA'jRaBflt) 

U20  ASSIGN  fPath  TO  I 

1630  SUBEUT 

1640  ' 

1650  ASSIGN  »Path  TO  Naiel 

1660  ENTER  ?Path  USING  "SA'iRaH*!)) 

1670  ASSIGN  fPath  TO  I 

1680  SUGEND 

1690  ' 

1700  ' 

1710  ! 

1720  SUB  Scope  setlRantOI) 

1730  OPTION  BASE  1 

1740  BEEP  1627. 60,. 3 

1750  INPUT  'ENTER  TINE  FOR  FIRST  DATA  POINT  lsec)',CliltO 

1740  INPUT  'ENTER  TINE  OF  EVENT  INITIATION  (secl'.EventO 

1770  INPUT  'ENTER  TINE  INCRENENT,  e.g.  lE-2  FOR  A  20  5  RUN  '.Tau 

1780  INPUT  "ENTER  VOLTAGE  RANGE  (V) '.Range 

1790  INPUT  "ENTER  DC  LEVEL  ZERO  FOR  CHANNEL  A  (VI',Dca 

1800  INPUT  'ENTER  DC  LEVEL  ZERO  FOR  CHANNEL  B  (VI'.Dcb 

1810  INPUT  'ENTER  BUCKING  VOLTAGE  FOR  CHANNEL  A  IV)',VaO 

1820  INPUT  'ENTER  BUCKING  VOLTAGE  FOR  CHANNEL  B  (V>",VbO 

1930  INPUT  'ENTER  SIGNAL  AHPLIFICATION  OF  ANALOG  COHPUTER',Aip 

1840    ! 

1850  Ra»«(ll='CURSR' 

1840  Raii«l2)=VflL«(Range/2000i           !  Vnon 

1870  ChO=(ABS!Chlt0l-ABS(Event0)lt2/Tautl 

1880  RaB«(3l=VALf(ChO) 

1890  Raii$(4l=VAL»(Taul 

1900  Rai<«(5l=VAL»(DcatlOOO)               '  Vczl 

1910  Rai(»(6)=VAL$IDcbll000)               '  Vcz2 

1920  Ra»«(7l=V«L»(Aiip) 

1930  Ra»»(ai=VAL«(VaO) 

1940  Raii«(9)=VAL«(ABS(VbO)) 

1950  SUBEND 

1960     ! 

1970    I 

1980     ! 


ii-..' 
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1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 

20  ?0 
2100 

:iio 

2120 

21  ;o 

2140 

2150 

2160 

2170 

21B0 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

22B0 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 


sue  Avg 
OPTldN  BASE  1 

con  /flrray/  Ra»(4096),fli512),fl(512),Ei512l,Zi512,4) 
con  /Datui/  [ii,Len.Ra.l'5,Ia,iau,I0,Va0,Vb0,Vczl,Vi::2,Vnor«,P« 
CON  /Stat;  ilip,«ode,Nav,Nblock,siqia 

:   THIS  SUBROUriNE  SEPARATES  THE  DATA  FOR  THE  STANDARD  RESISTOR 
I   PLUS  THE  WIRE  {AdI}  AND  THE  HIRE  {BID)  AND  CALCULATES 
'   THEIR  RUNNIN6  AVERAGES 

IF  l1oi)e=2  THEN  Ch5tart=Ra»(IO)-l 

J5tart=Ch5tartH 

FOR  1=1  TO  Nblotk 

Suii=0 

SuaY=0 

Jerd=Jstart«2IHav-l 

FOR  J=Jstart  TO  Jetid-l  STEP  2 
Sum^Sutxi-Raii  IJ I  tVnor*/A>p 
Suny^SuiY+RaKlJ^I 1 (Vnors/ftio 

NEXT  J 

A(I)=Suis/Nav-Vczl/A<ipi-VaO 

B(I)=Sui>y/Nav-Vcz2/AiiptVbO 

Jstart^Jendtl 
NEKT  I 
5UBEHD 


CONVERT  TO  TRUE  VOLTAGE 
AND  AVERAGE 


SUBTRACT  ZERO  LEVEL  AND 
ADD  BUCKING  VOLTAGES 


SUS  Res is 
OPTION  BASE  1 

CON  /Array/  RaH(4096),A(5121,Bi512),EI512),i(512,4) 
CON  /Oatuii/  Di,Len,Ra,R5,Ta,Tau,T0,Va0,Vb0,Vc:l,Vc22,Vnori,P« 
CON  /Stat/  Asp, Node, Nav,Nblock,Siq«a 

THIS  SUBROUTINE  CALCULATES  THE  INITIAL  RESISTANCE  OF  THE  HIRE 


IF  Node<>2  THEN  SOTO  2381 
Ra=-VbO/VaOIRs-. 00277 
3UPEX1T 


Suix=0 
Suix2=0 
Suiy=0 
Su«y2=0 
Suiixy=0 
iz=10 
FOR  1=3 
JI=(I- 


FOR  INITIAL  FLUK  ON  WIRE 


NO  INITIAL  FLUX,  USE  LEAST  S8UARES  FIT 


TO  12 
.5)INav 

R=-B(I)tRs/AiI)-.0O277 

Tiie=JI"1.25 

Suix=Suflx>Tiie 

Su«x2=Suix2Hiiie"2 

Spty=Suily+R 

Siny2=Suiy2tR"-2 

Suixy=Siilixy+TiiietR 
NEXT  I 

Betal  =  (Suisy-(Sui!lSu«yl/!:)/(Si]ix2-(Su«x'-2)/Z!) 
BetaO=(Suiy/7zl-BetalllSuix/Z!) 
S=Sii«y2-(Suiy'2l/Zz 
T=Suix2-ISu«x"21/Zz 
U=l./(Zz-2.l 
V=Beta|-2 
N=UtlS-VIT) 

Sigta=SeR(HISuiix2/Zz/TI 
Ra=fletaO 


TINE 

ENPIRICAL  ADJUSTHENT  FOR  BETTER  FIT 


VARIANCE  ABOUT  REGRESSION 
EST  OF  STD  DEV  OF  INTERCEPT 


BEEP  2197. 26, ,3 

PRINT  ■  ■,"RA  =  'iRa;"  AT  'iTO!"  C 

Neii=Ra 

INPUT  "Enter  no»  value  for  RA  or  accept  default", New 

IF  Ne»=Ra  THEN  2720 

PRINT  'III  NEN  VALUE  ltr,"Rfl  =  •;Ne»;'AT  ";T0;"  C 

Ra=Neit 
SUBEND 
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2730  ! 

2740  ! 

2750  ' 

2740  SUB  Props 

2770  OPTION  BASE  1 

2780  COH  /ftrray/  Ra«f4096),fli5l2),B(512),E(512),Z(512,4) 

2790  CDN  /Const/'  Alph,Beta,flbeg,l(cnil 

2800  con  /Datum/  Di,Len,Ra,R5,Ia,rau,T0,Va0,Vb0,Vr!l,Vr22,Vnor«,P« 

2810  CDN  /Stat/  fl«p,l1odE,Nav,Nbiock,Sig«a  ' 

2920  ! 

2830  '   THIS  SUBROUTINE  CALCULftTES  THE  PROPERTIES  OF  THE  HIRE: 

2840  : 

2850  !    i(.],I)  =  TIHE 

2860  !    Z(J,2)  =  SUPERFICIAL  HEAT  FLUK 

2870  !    n.J,31  =  THETA  [  T(iiire)  -  Tlanbientl  ] 

2B80  !    Z(J,4)  =  NGBIFIED  NUS5ELT  NUMBER 

2870  ' 

2900  CDn5t=fia/(I.tAlphlT0*BetalT0'2)        !     Rlitirel  AT  T=0  C 

2910  IF  l1ode=l  THEN  2990 

2920  Neo=INT({Raitl3)-RaiillO))/Na¥/21 

2930  PRINT  "  ■/HIRE  PROPERTY  i-/R(»ire)  I  AT  T=0  C  I  =  'iConst 

2940  Ne»=Const 

2950  INPUT  "Enter  new  value  for  Rliiirel  or  accept  default", Ne« 

2940  IF  Neii=CDn5t  THEN  2990 

2970  PRINT  'III  NEN  VALUE  m","R(iiire)  [  AT  T=0  C  I  =  "iNen 

2980  CDn5t=Neit 

2990  Bb=AIqh/Beta 

3000  Ba=-.5IBb 

3010  B2=Bb"2 

3020  FOR  1=1  TO  Nblock 

3030  ZII,ll=TaulNavl(I-Neg-.5) 

3040  Z(I.2)=AiI)tA(I)l(-B(I)/A(I)-.00277/R5)/(R5»3.14159IDilLen)-Qbea 

3050  R=(-B{IIIIRs/A(n-. 00277 

3040  F=B2t!4./Beta)l(R/Con5t-l.l 

3070  Z(I,3)=Ba-.5IS8R(FI-Ta 

3080  Z(I,4)=Z!I,2)IOi/acndlZ(I,31) 

3090  ' 

3100  ErI=lVnori"2)/iA(II'2) 

3110  Er2=(Vnnrii'2)/(B(I)"2) 

3120  Er3=(Siq«a'2)/IRa''2l 

3130  EiI)=IErltEr2tEr3l/((l-(Rj/R)l*2l 

:i40  NEXT  I 

3150  SUBEND 

3140 

3170 

3180 

3190  SUB  0utp(N«, Title!) 

3200  OPTION  BASE  1 

3210  CON  /Array/  Ra«(4096) ,A(512),Bi5121,E(512) ,Z(512,4) 

3220  COH  /Datui/  Bi,Len,Ra,Rs,Ta,Tau,T0,Va0,Vb0,Vczl,Vcz2,Vnori,P« 

3230  COH  /Stat/  Anp. Node, Nay. Nblock, Siqia 

3240  DIM  Addf[161,Hel»[aOI,He2l[BO!,He5s[80I,Inb«t71,PrB«[181,Pl«[80I,Rl$I801 

3250  DIN  St«[80],SI([801,IlJI801,KlJ[80I,H2«t80],Outl(3ll80],a(4)[loi 

3260 

32/0  Npg=!NI(Nblock/47)             !  INFORHATION  ABOUT  FILE  SIZE 

3280  Ni=Nblock-Npgl47 

3290  IF  Nx>42  IHEH  NbIork=Nblock-iNs-42) 

3300  IF  N!!=0  IHEN  Npg=l(pq-1 

3310  Npg=Npgt| 

3320  5ize=I9INpo 

3330  ALLOCATE  6iNblork,2!     '  "BUICK  i  DIRTV"  PLOT 

3340  Abs=3 

W3.  ?r''i^.  ,  ,„.„  ..      '   "  "EAT  FLUX  V  THETA   if  PRESSURE  =  const 

3o40  IF  ;ode=2  THEN  Abs=l    '   -  THETA  v  TIHE      if  HEAT  FLUX  =  const 

3o70  IF  l(ode=2  THEN  0rd=3 

3380  FOR  1=1  10  Nblock 

3390  G(l,l)=Z(l,Ab5! 

3400  G(I,2(=ZiI,0rd) 

3410  NEXT  I 

3420  ' 

3430  BEEP  1427. 40,. 3 

3440  PRINT  CHR«(l2i 

3450  CALL  Scr  plot(6fn) 

3440  DEALLOCATE  6(li 
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:!47(')  ' 

3480  INPUT  "Do  vou  msh  to  have  a  HARD  COPY  of  the  nuitiers  ?  I  y/N  ).',fln5« 

3490  IF  AnsJO"!"  THEN  4060 

3500  PRINT  TitieJ,'  ","  ■/SOURCE  DATA  FILE:  •■,!» 

3510  PRINT 

3520  PRINT  ■tilt     PLEASE  INSERT  A  DISC  INTO  DISC  DRIVE  1' 

3530  PRINT 

3540  PRINT  -AN  ASCII  FILE  KILL  BE  CREATED  ON  THE  DISC.  THIS  FILE  CAN  THEN^ 

3550  PRINT  "BE  TRANSFERED  TO  THE  HAIH  COHPUTER  FOR  PRINT  OUT." 

3560  INPUT  "PLEASE  ENTER  A  UNIQUE  NAME  FOR  THE  ASCII  FILE.',A5« 

3570  Pre«=' 

3580  Add«=" 

3590  Int)*=^ 

3600  St«='           tutttttttttlttttltltttlttlttlltlttttttlttttlttttf 

3610  HeH=^            TINE    HEAT  FLUX    THETA      NU55ELT" 

3620  He2«='             is)      (H/inli)      fC)       NUMBER" 

3630  Hp3$="                     page  " 

3640  Ra«=VAL«iRa) 

3650  Ra$=HaUl,7] 

3660  PlMAifdH'PRESSURE     =  "SPU"  kPa" 

3670  Rl«=AddH'INITIAL  RESIS.  =  ■4Ra»4^  OHNS  AT  ■JVALKTOIl^  C" 

3680  5iq«a«=VAL»(Siqiial 

3690  Sig«a»=5iq«anl,7] 

3700  SIMddtV'SIGUA       =  "4SigH« 

3710  Tl«=flddirAHBIENT  TEKP.  =  "iVALMTali"  C" 

3720  KI»=AddH"HIRE  LENGTH   =  "iVALXLenll^  «■ 

3730  K2«=Add»S"HIRE  DIANETER  =  'WALJIDi)!.'  N" 

3740  Np9=l 

3750  Lin=l 

3760  NABS  ST0RA6E  IS  ■;HP82901, 700,1" 

3770  CREATE  ASCII  AsJ.Siie 

3730  ASSIGN  fPath  TO  kt 

3790  OUTPUT  fPath:Add«iTitIe»,Inh«,NI».K2$ 

3800  OUTPUT  ?Path!Tl«,Pl«,Rl« 

3810  OUTPUT  SPath;Sl$,Irtb»,St«,lnb«,Hel$ 

3820  OUTPUT  »PathjHe2i,Inbi,St*,Inb* 

3830  FOR  1=1  TO  Nblocl;  STEP  3 

3840  CALL  Dec(I,Z«(ll] 

3850  0ut«(ll=Pre«4Z«(nilnb»it$(2)»Inb»SZ$(3)!!lnb«U»l4) 

3860  IF  l+l':=Nbloci<  THEN  CALL  Dec(Itl,Z«(l)) 

3670  Gut«(2l=Pre(!,Z«ll)iInbHZJ(2IHnb«iZ«(3)<(Inb«Z«(41 

3880  IF  It2<=NbIock  THEN  CALL  Dec(It2,Z«lll) 

3890  0ut$(3l=PreHZ»lllJInb«iZll2)iInb«i,Z«(3)!iInbnZI(4l 

3900  IF  (ItlDNblock  THEN  0utJ(2l=AddJ 

3910  IF  (lt2)>Nblock  THEN  Oijtll3)=Add» 

3920  OUTPUT  »Path;0ut«ll),0ut$(2l,0ut$(3l 

3930  Lin=Lintl 

3940  IF  Npg=l  AND  Lin=15  THEN  SOTO  3960 

3950  IF  Lin=17  THEN 

3960  OUTPUT  PPath:Inb»,Inb«,Inb«,Inb«,Inb«.Inb$,Inb(,Inb« 

3970  Tl«=Add«4Title«iHe3«VAL«INpoH) 

3980  OUTPUT  ePathiTl«,Inb$,St«,Inbt 

3990  OUTPUT  fPathiHel«.He2«,Inb«,Stl,Inb»,Inb» 

4000  Lin=l 

4010  Npg=Npgtl 

4020  END  IF 

4030  NEXT  I 

4040  flASS  STORAGE  IS  ■:HP82901,700,0^ 

4050  BEEP  1302. 06,, 3 

4060  I 

4070  INPUT  "Do  you  wish  to  have  a  PLOT  of  the  nuibers  (Y/NI  ?",Ans« 

4080  IF  Ans«<>"Y"  THEN  4610 

4090  PRINT  CHR«(I2) 

4100  PRINT  Title*, ■  ","  ", "SOURCE  DATA  FILE:  ";N« 

4110  PRINT 

4120  PRINT  'till     PLEASE  INSERT  A  DISC  INTO  DISC  DRIVE  1" 

4130  PRINT 

4140  PRINT  "BOAT  FILES  NILL  BE  CREATED  ON  THE  DISC,  THESE  FILES  CAN  THEN" 

4150  PRINT  "BE  PLOTTED  USING  THE  GENERAL  PLOTTING  PROGRAM" 

4160  PRINT 
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tl70 
4180 

4H0 

4200 
4210 
4220 
4230 
4240 
4250 
4240 
4270 
4280 
4270 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4360 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4B70 


PRINT  "IHE  TYPE  DF  PLOTS  POSSIBLE  ARE:" 

PRINT 

PRINT  ■  1  =  NUSSELT  KUHBER  vs  TIKE" 

PRINT  "  2  =  HEfll  FLUJ  v5  IHETA" 

PRINT  ■  3  =  BOTH  PLOTS  1  «ND  2   (plot  1  retoriipd  first)" 

PRINT  "  4  =  HEM  FLUK  vs  TIHE" 

PRINT  ■  5  =  THEIfl  vs  TINE" 

PRINT  "  6  =  BOTH  PLOTS  4  «N0  5 


INPUT  "CHOOSE  THE  TYPE  OF  PLOT  YOU  HISH. 

IF  Hode=l  THEN  CALL  Hani 
Foriat»="(1Z.BDDE,r 
ALLOCATE  BiNbli)tk,2) 
(lASi  STORAGE  IS  ":HP82901,700, 1' 
IF  9()=!  OR  Bq=3  THEN  4340 
IF  Sq=2  THEN  4370 
IF  Qq=4  OH  8q=6  THEN  4400 
IF  6q=5  THEN  4430 

Abs=l 

OrdM 
SOTO  4460 

Abs=3 

0rd=2 
6010  4460 

AbsM 

Ord=2 
60TO  4460 

Abs=l 

0rd=3 


INPUT  "PLEASE  ENTER  A  UNIQUE  NflHE  FOR  THE  BDAT  FILE, 
FOR  0=1  10  Nblock 


iplot  4  recorded  first)" 
,Qq 


,Bd$ 


G(0,l)=ZIO,Abs) 

SI0,2)=Z(0,0rd) 
NEXT  6 

CREATE  BDAT  Bd«, l2)Nblock).9 
ASSIGN  PPath  TO  Bd« 
OUTPUT  fPath  USING  ForiatljGltl 
ASSIGN  fPath  TO  I 
PRINT 

PRINT  "BDAT  file:  "iBdJiNblockj"  points', "Data  foriat:  HZ.DDDE" 
IF  Bq=3  THEN  4370 
IF  Sq=6  THEN  4430 
DEALLOCATE  Gdl 

MASS  STORAGE  IS  ":HP82901, 700,0" 
BEEP  5126.94,.! 
PRINT  "    END" 
5UBEND 


SUB  Decd.Ztit)) 
OPTION  BASE  1 
COM  /Array/  RaM(4096),A(512) ,B(512) ,E(5i2) ,ZI512,41 

THIS  SUBROUTINE  ENSURES  THAT  ALL  VALUES  HAVE  THE 
CORRECT  FIELD  LENGTH  FOR  OUTPUT 

Lonq=5 
Dec«="0" 
Perl="." 
FOR  1=1  TO  4 

Zr=Z(l,J) 

IF  Z:(0  THEN  Zz«=VALt(DR0UND(FNRound(Zz.3),5)) 

IF  Zz<0  THEN  4880 

Z2i=VAL»(DR0UND(Zz,5)) 

IF  Z:<1  THEN  4930 

6010  4980 

Z)(J)=Z:» 

Lonq=6 
NEH  3 
GOTO  5070 
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4880 

4B90 

4900 

4910 

4920 

4930 

4940 

4950 

4940 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

508(1 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230. 

KM; 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5120 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

550(1 

5510 

5520 

5530 

5540 

5550 

5540 

5570 

5580 


FOR  Zi<0 
IF  LEN(2!«l>Long  THEK  i:«=Zzltl.Loiig] 
IF  LENIZzlKLong  THEN  Zz»=Z!$l.Dec* 
IF  LENIZzJIOLonq  THEM  4900 
SOTO  4840 

FOR  OiZZfl 
IF  LEN(Z2«)>Long  THEN  Zz«=Z:»Il,Ung-n 
IF  LENIZ2»KLi)ng  THEN  Z:«=Dec«4Z:« 
IF  LENlZzSIOLong  THEN  5040 
GOTO  4840 

FOR  ZZ>1 
P»="N" 
FOR  F=LEN(Zz«l  TO  1  STEP  -1 

IF  Z2»IF,F]=',-  THEN  P$="V 
NEIII  F 

IF  P«='N"  AND  LENIZztXLong  THEN  Zz«=Zz«iPer« 
IF  Zz<l.E*5  AND  LEMIZzfKLong  THEN  Zz»=Zz»!iDec« 
IF  Zz<l.Et5  m  LEH(Zz«)<>Long  THEN  5040 
GOTO  4840 
SUSEND 


DEF  FNRoundlZ.Np) 
Nu=INTIZI10'Npt.5) 
Nu=Nu/(10"Np) 
RETURN  Nu 

FNEND 


ROUNDING  FOR  Z<0 


SUB  Nani 
OPTION  BASE  1 

con  /Array/  Ra»f40961 ,S(512),B(512),E(512),Z(512,4I 
ton  /Stat/  Aip,node,Nav,Nblocl:,Sigia 

THIS  SUBROUTINE  EVALUATES  THE  ERROR  IERU  CALCULATED 
IN  SUBROUTINE  PROPS  IN  ORDER  TO  EL111INATE  THOSE  FIRST 
VALUES  THAT  HAVE  5!  OR  HOPE  ERROR 

Ht=0 

FOR  1  =  1  TO  Kblocl! 

Ht=I 

IF  (EIIX.05  AND  Z(I,3)>0.I  THEN  GOTO  5330 
NEIT  1 

Nblock=Nblock-Nt 
IF  Nblock=0  THEN  GOTO  5410 
FOR  .J  =  l  TD  Nblocfc 

ZIJ,U=ZIJ»nt,l) 

Z(J,2)=Z(J*nt,2) 

ZIJ,3)=ZIJtHt,3l 

ZIJ,4l=Z(JtMt,4) 
NEJT  J 
SUBEND 


SUB  Graf  setlTscalel 
con  /Sht/  Aip,llDde,Nav,Nblpck,Sigaa 

READ  X$,y«,)(!il,l(x2,¥yl.yv2 

IF  node=2  THEN  READ  U,VJ,hl,l(x2,Vyl,¥y2 

IF  I1p(le=2  THEN  X»2=T5cale 

BINIT 

GRAPHICS  OFF 

VIEKPORT  0,130,15,100 

PEN  1 

LDIR  PI/2 

LORS  6 

mVE  0.57 

LABEL  '/« 
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5590  LDIR  0 

5o00  L0R6  4 

5410  HOVE  65,15 

5620  LABEL  U 

5630  CSIZE  3. 7, .5 

5640  IWVE  15,15 

5650  LhBEL  Jil 

5660  MOVE  127.15 

5670  LABEL  X«2 

5680  LOPS  2 

5690  ilOVE  0,20 

5700  LABEL  Vyl 

5710  HOVE  0,98 

5720  LABEL  Yy2 

5730  IF  ((ode=2  THEN  5780 

5740  Xx1=LGTIXk1) 

5750  l(ii2=L5I(ll)(2) 

5760  Vyl=L6T(VyI) 

5770  Yy2=LGTiyy2) 

5780  VIEKPDRT  15,130,20.100 

5790  KINDOU  l(xl,Kx2, Wl,Vy2 

5800  FRABE 

5810  IF  node=l  THEN  5840 

5820  HOVE  0,yyl 

5830  DRAH  0,Vy2 

5840  ' 

5850  DATA  "THETA-.-HEAT  FLUr',10,100,lE4,lE6 

5860  DATA  '  TINE',"  THETAV-2,20,30,70 

5870  SUBEND 

5BB0 

5890 

5900 

5910  SUB  Scr  plot  (SO)  I 

5920  OPTIOR  BASE  1 

5930  con  /Stat/  A«p,node,Nav,Nbl!jck,3iqia 

5940  Er$="N' 

5950  l5tart=l 

5960  IF  l1ode=2  THEN  6030 

5970  FOR  1=1  TO  Nblock 

5980  IF  8(1,1X0  OR  811,2X0  THEN  l5tart=IH 

5990  IF  htart=I*l  THEN  6020 

6000  S(1,11=L6T(G(1,1II 

6010  6(1,21=LGT(S(1,2)) 

6020  NEIT  I 

6030  LINE  TYPE  1 

6040  PEN  1 

4050  ALPHA  OFF 

6060  GRAPHICS  ON 

NOVE  5(Istart,ll,6(l5tart,2) 

FOR  I=htart  TO  Nblock 

6090  DRAN  G!1,1I,6(I,2) 

6100  NEXT  1 

6110  HOVE  0,0 

6120  IF  Er«='Y"  THEN  6190 

6130  IF  «ode=l  THEN  DISP  "     III  LOG-LOG  SCALING  I  I  I" 

6140  NAIT  5 

6150  INPUT  ■     ERASE  THIS  PLOT  (Y/N)  '',Er» 

6160  IF  ErSO'Y'  THEN  6190 

6170  PEN  -I 

6180  GOTO  6070 

6190  GRAPHICS  OFF 

6200  ALPHA  ON 

6210  SUBEND 


6070 
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APPENDIX  B 
Program  PRESSURE:  Pressure  Data  Analysis 

This  program  was  written  in  Hewlett-Packard  BASIC,  Version  2.0,  and 
run  on  the  HP  9816  microcomputer.   Experimental  parameters  are  Included 
in  the  voltage  input  data  set,  similar  to  the  pressure  transient  option 
of  program  BOIL.   Manipulation  of  data  provided  by  the  digital 
oscilloscope  is  as  described  in  Appendix  A;  however,  since  the  entire 
capacity  of  the  oscilloscope  is  used  to  record  the  pressure  signal, 
separation  of  data  by  input  channel  is  not  necessary.   The  program 
structure  is  identical  to  that  of  BOIL,  and  the  output  is  arranged  such 
that  instantaneous  pressures  are  reported  at  the  same  experimental  times 
as  the  output  from  BOIL. 
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ItlttttttlOtUtttttttttttt  PPESSUFiE  tttntltlttlttttttUttllt 

II 

PflOGRflK  CONVERTS  VOLTflSE  DATA  ilEflSURED  FROM  fi  PRESSURE 
TRANSDUCER  (VIA  A  CHARGE  ANPLIFIER)  TU  A  PRESSURE  HISTORY 
OF  THE  PRESSURE  VESSEL.  VOLTAGE  DATA  IS  ORIGINALLV  INPUT 
FROn  THE  NICOLET  DIGITAL  OSCILLOSCOPE,  BUT  THESE  DATA  ARE 
ALSO  STORED  ON  A  DISC  FOR  FUTURE  REFERENCE. 
-->NOIE<:--   D-SCOPE  IS  CURSOR  TRIGGERED. 


D.SCHNIDT 


m       SEP  1984 


update:  'quick  i  iirty'  plot  (D.SCHNIBT)  KSU   JAN  1985 


itDiitttttitiiitiuunntititiutitttuttuttitttiutuittuttt 

1 

OPTION  BASE  1 

con  /Stst/  ChO, Drift, Nav.Nblock,Pperv,Ta.Iau,Vm)r« 

CON  ,'Array/  A(1025),RaH(4096l 

DIN  Title«[80I 

PRINT  CHR$(12) 

INPUT  "RESET  GRAPHICS  lY/N)  ?', Reset* 

INPUT  -NHAT  IS  THE  RUN  NUNBER  ?', Title* 

Title«="RUN  NUNBER  ■iTitle* 

PRINT  Title* 

HA IN  PROGRAN 

CALL  Inpt(N$) 
T5i:ale=Taul4000 

IF  Re5et«='V"  THEN  CALL  Graf  5et(RaN(6),Tscale) 
PRINT  Title*, ■  ",■  ", "SOURCE  DATA  FILE;  •■,!» 
PRINT 
PRINT  "AVERAGING  DATA  >  "iNav;"  POINT  AVERAGING" 

CALL  Avo 
PRINT  "CALCULAIING  AHP  DRIFT" 

CALL  fltp  drift 
PRINT  "CALCULATING'PRESSURE  HISTORY" 

CALL  Pressure 

CALL  Oiitp(NS,TitIe»l 
INPUT  "RECALCULATE  NITH  SANE  SOURCE  DATA  (Y/N)  ?",Aiis* 
IF  An5«>"Y"  THEN  400 
PRINT  CHR«(12) 
GOTO  270 
END 


SUB  Inpt(N*l 
OPTION  BASE  1 

CON  /Stat/  ChO,Drift,Nav,Nblock,Pperv,ra,Tau,Vnor« 
CON  /Array/  A(1025),Raiil4096) 
DIN  Raii*l4096iI51 

r 

INPUT  "DATA  SOURCE  SCOPE  =  1,  BOAT  FILE  =  2",Tt 

CALL  Reader (Tt,N»,Ra«*(l) I 
FOR  1=2  TO  4096 

Ra»(Il=VflL(Ra»*(II) 
NEXT  I 

Vnori=Ra»(2) 
ChO=Rai((3l 
Iau=RaK(4) 
Ta=Ra»(51 
BEEP  2197. 26, .3 
Nav=8 

INPUT  "ENTER  [4,8,161  POINT  RUNNING  AVERAGE  {DEFAULT'S!  ",Nav 
Nch=4096-ChO-Na¥/2 
Nblock=INTINch/Nav) 
PRINT  CHR*I12) 

PRINT  Rai(*(ll,Ra»l2),RaH(3l,RaK(4),Ra«i(5),Raii(6),Rai((7) 
PRINT 
SUBEND 
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SUB  Rea[ler(Tt,Na«e«,Ra«S(tl) 
OPTION  BfiSE  1 
PRINT 

PRINT  "PLEASE  INSERT  DftTA  DISC  INTO  DISC  DRIVE  0" 
INPUT  -PLEASE  ENTER  fl  UNIQUE  NfillE  FOR  THE  BD«T  FILE.",Na«e» 
DISP  "HORtING,  PLEhSE  WIT." 
NflSS  STORAGE  15  " :HP82901, 700,0" 
IF  Tt=2  THEN  1000 

ASS16N  sScQpe  TO  9 
CONTROL  9,3i9600 
ASSIGN  fScopeiFORNAT  ON 
CONTROL  9,5!3 
CONTROL  9, 4; 2 
OUTPUT  ?5cope;CHR»(l); 
OUTPUT  »Scope;CHR>(69l;CHR$(48)|CHR«(68l;CHR«(49);CHRI(68);CHR«(48): 
OUTPUT  »Scope;CHR«(79);CHR$(52)iCHR»(43);CHR>(57)iCHR$l54) 
OUTPUT  «5copeiCHR»(2); 
ENTER  fScppe  USING  "5fl,ll,X";RaH»(l) 
CONTROL  9,5i0 
ASSIGN  "Scope  TO  I 

CALL  Scope  set (Raii«(l)l 

DISP  "NORKINS,  Please  nait." 

CREATE  FDAT  Na«e«,4096,5 
ASSIGN  ePath  TO  Naiet 
OUTPUT  ePalh  USING  "SA'iRaBXIl 
ASSIGN  »Path  TO  I 
SUBEXII 

ASSIGN  fPath  TO  Naie« 
ENTER  ePath  USING  "SA'iRaiitd) 
ASSIGN  ePath  TO  I 
SUBEND 


SUB  Scope  setlRantdll 
OPTION  8ASE  1 
BEEP  1627. 60,. 3 

INPUT  "ENTER  TINE  FOR  FIRST  DATA  POINT  l5ecl",ChltO 
INPUT  "ENTER  TIME  OF  EVENT  INITIATION  (sec)",EventO 
INPUT  "ENTER  TINE  INCRENENT,  e.g.  5E-3  FOR  20  SEC  RUN",Tau 
INPUT  "ENTER  VOLTAGE  RANGE  (V)",Ranqe 
INPUT  "ENTER  INITIAL  PRESSURE  (psigCPo 
INPUT  'ENTER  FINAL  PRESSURE  lp5ig)',Pf 
INPUT  "ENTER  ANBIENT  TEHPERATURE  (C)",Ta 

Raiillll="PRESS"  I  For  error  if  read  bv  "BOIL" 

Raii«(2)=VAL$(Ranqe/20001  !  Vnom 

Ra\(»l3)=VAL$((ABS(ChU01-ABS(EventO))/Tau)   !  ChO 
Raii«(4)=VAL«!Tau) 
Ra»«(5)=VALS(Ta) 
Rai(>(61=VAL$(Pol 
Rfl»«(71=VAL»IP() 
SUBEND 


SUB  Avg 
OPTION  BASE  1 

con  /Stat/  ChO, Drift, Nav,Nblock,Pperv,Ta,Tau,Vnori 
con  /Array/  A(1025),Rai((4096l 

THIS  SUBROUTINE  CALCULATES  THE  RUNNING  AVERAGES  OF  THE  DATA 

A(li=lNT(Nav/2)-l 
jBtart^ChO-Adj 
FOR  1=1  TO  Nblock 

Suii=0 

Jend'JstarttNav-l 

FOR  J=j5tart  TO  Jend 
Suii=Su«tRa»(J)tVnori« 

NEU  1 

A(I)=Suii/Nav 

J5tart=jendt| 
NEXT  I 
SUBEND 


CONVERT  TO  TRUE  VOLTAGE 
AND  AVERAGE 
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H80  ! 

149f!  ' 

1500  ! 

1510  SUB  Aip  drift 

1520  OPTIOB  BASE  1 

1530  CON  ,'Stat/  ChO, Drift. Nav,Nbiiicl;,Pperv, la. Tau.Vncri 

1540  con  /ftrrav/  (1(1025), Raiii4096) 

1550  ' 

1540  '   THIS  SUBROUTINE  CflLCULflTES  THE  DRIFT  IN  THE  CHftRSE  m  SIGNAL 

1570  '   IT  IS  ASSUMED  THAT  DEPRESSURHflllON  IS  COHPLETE  BV  THE  HUE 

1580  '   OF  INTEREST  FOR  THIS  CALCULATION 

159(1  ' 

1600  Sii«]i=0 

1610  5u««2=0 

1620  Sij«y=0 

1630  Siiixy=0 

1640  Z:=H 

1650  FOR  I  =  (Nblocli-15)  TO  (Nblock-2) 

1660  T=IautNavl(I-I)                  I  TIKE 

1670  SuiX=Sii«xtT 

16B0  Su««2=Suix2tT'2 

1690  5uiy=Su«ytA(I) 

1700  Suiixy=5iiii«vtAII)tT 

1710  NEXT  I 

1720  Drift=DR0UND((Su«xy-(SuiixlSii»yl/Zz)/(Su«x2-(Sij«x'2)/Z:),6) 

1730  SUBEND 

1740 

1750 

1760 

1770  SUB  Pressure 

1780  OPTION  BASE  1 

1790  con  /Stat/  ChO,Drift,Nav,NbIock,Pperv,Ta.TaD.Vnorii 

1800  con  /Array/  A(1025l,Raii(4096) 

1810  Din  Test  (10251 

1G20  ! 

1830  1=0 

1840  A.in=100 

1850  FOR  1=1  TO  Nblock                  I  ADJUST  VOLTAGE  FOP  DRIFT 

1860  Vdrift=TaulNa¥l(I-l)»Drift 

1870  Te5t(I)=A(I)-Vdrift 

1880  IF  TestdOAnin  THEN  1910 

1890  A«in=Test(ll 

1900  Iiin=I 

1910  NEST  1 

1920  Pper¥=DR0Umi(Rsi((6)-Raw(7)l/(A(ll-Aiin),6)    >   psi/valt 

1930  ' 

1940  IF  m   THEN  BEEP  2197. 26, .3 

1950  PRINT-  ", "DRIFT  =  'jDrift;'  Volts  per  second','  ■,'PSI/VOLT  =  "jPperv 

1960  Nei(=Dri(t 

1970  INPUT  "Enter  rei<  value  for  DRIFT  or  accept  default", New 

1980  IF  Neii=Drift  THEN  2030 

1990  Drift'NsK 

2000  1(=1 

2010  DISP  "RECALCULATING  — >  PSI/VOLT' 

2020  GOTO  1840 

2030  FOR  1=1  TO  Nblocl! 

2040  A(II=Ra»(4)-(Iest(l)-Test(I))IPpery         '  CONVERT  TO  PRESSURE 

2050  IF  Dliir  THEN  A(Il=A(Iiin) 

2060  NEIT  I 

2070  SUBEND 

2080 

2090 

2100 

2110  SUF  Outp(N»,TitIe«) 

2120  OPTION  BASE  1 

2130  CON  /Stat/  ChO, Drift, Nav, Nblock, Pperv,Ta,Tau,Vnori 

2140  con  /Arrav/  A(1025),llaN(40?6l 

2150  Din  Add«[l6],HelH8()],He2«I80I,He3«[B0],Inb«(5],St$[8Ol 

2160  DIN  Tl«tBOI,T2»IB01,Out«{3)[801,Z«(4![101 

2170  I 

2180  Istart=I6/Na«H                   !  FILE  SIZE  INFORMATION 

2190  Step=(Istart-l)l2 

2200  Si:e=INT(Nblock/Stept.5i 

2210  ALLOCATE  G(5i:e,2) 

2220  6(1,1)=0 

2230  G(l,2)=Raiil6) 
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2240  J=l 

2250  FOR  I=Istart  ID  Nblock  STEP  Step 

2260  J=J*1 

2270  IF  .]>Si:5  THEN  2300 

2280  G(J,11=TjiiINjvI(I-1)               !  IIHE 

2290  6U,2)=ftlI)                    !  PRESSURE 

2300  NEST  I 

2310  ' 

2320  PRINT 

2330  PRINT  ■  ■/ItUtttlltltltUttUDIItttttllilttltltttlinitlltttf 

2340  PRINT"  ",■     the  "QUICK  !(  BIRIV""  screen  plot  will  be" 

2350  PRINT  ■  ","   F516  v  TIKE  regardless  of  output  units  choice" 

2340  INPUT  "ENTER  PRESSURE  UNITS  FOR  OUTPUT  !  f.  =  kPa,  P  =  p5ig",Unit» 

2370  BEEP  1427.60. .3 

2380  PRINT  CHR$(121                    !  "QUICK  I  DIRTV  PLOT 

2390  CALL  Scr  plot (S(l),Si!el             '  PSI6  v  TINE 

2400  I 

2410  IF  Unit«="P"  THEN  2470 

2420  DISP  "CONVERTINS   psig >  kPa  ■ 

2430  Pperv=Pperv/14.496«I01.325 

2440  FOR  1=1  TO  Siie                  !  CONVERT  TO  kPa  (absolute) 

2450  6(1, 2)=i6(l,2)/14. 694*111101. 325 

2440  NEJT  I 

2470  INPUT  "Do  you  »ish  to  have  a  HMD  COPY  of  the  nuibers  ?  (Y/N)",flns» 

2480  IF  «ns«<)"i"  THEN  3000 

2490  PRINT  Title*,"  ","  ", "SOURCE  DATA  FILE:  "iN« 

2500  PRINT 

2510  PRINT  "till     PLEASE  INSERT  A  DISC  INTO  DISC  DRIVE  1" 

2520  PRINT 

2530  PRINT  "AN  ASCII  FILE  KILL  BE  CREATED  ON  THE  DISC.  THIS  FILE  CAN  THEN" 

2540  PRINT  "BE  TRANSFERED  TO  THE  NAIN  COHPUTER  FOR  PRINT  OUT" 

2550  INPUT  "PLEASE  ENTER  A  UNIQUE  NAHE  FOR  THE  ASCII  FILE  ".flsl 

2540  Add«=" 

2570  Inb«=" 

2580  St«="         IttttltlttttUtlttttllt       IttllttttltlltlttDtItt" 

2590  Hel«="           TINE    PRESSURE          TINE    PRESSURE" 

2400  He2$="            (s)     (psig)            (s)     (psig)" 

2610  He3«="                   PAGE  " 

2620  Tl«=Add«!,"AMBlENT  lEHP.  =  "tVALJITa)*"  C" 

2630  T2«=HddH"PSl  PER  VOLT   =  "lVAL«lDR0UND(Pperv,6)) 

2640  IF  Unit«="P"  THEN  2670 

2450  He2«="            (s)     (kPA)             (s)     (tPa)" 

2660  T2«=Addllc"kPa  PER  VOLT   =  "!iVAL«(DR0UNDiPperv,4)) 

2670  l5tart=l 

2480  lendMS 

2490  Npg=l 

2700  NASS  STORAGE  IS  ■:HPB2901, 700,1" 

2710  CREATE  ASCII  A5«,60 

2720  ASSIGN  »Path  TO  kt 

2730  OUTPUT  fPath;Add«lcTitIe«,Inb«,TU,T2».Inb« 

2740  OUTPUT  »Path;St«,Inb«,Hel«,He2« 

2750  OUTPUT  sPathjInbl.StJ.InbJ 

2740  FDR  I=l5tart  TO  lend  STEP  3 

2770  CALL  Deed, 6(1), Ztdl, Size) 

2780  0ut«(l)=Add«Z»(I)4!nh«ui(2)!rAdd$lZ»(3)tInb«lZ»(4l 

2790  IF  Itl<=Size  THEN  CALL  Dec(I*l,6(tl,i«(ll,SizB) 

2S0O  Out«(2)=AddHZ$(l)lInb«!!Z«(2)iAdd$!(7l(3)HnbJliUI4) 

2810  IF  I*2<=Size  THEN  CALL  Dec(i*2, Gil), Z$(l), Size) 

2820  0utJ(3)=Add$i.Z»(l)llnb*iZ«(2)!(Add«!iZI(3)!'inb«l!Z$(4) 

2B30  IF  I+DSize  THEN  Out«(2l=Add$ 

2840  IF  l*2>Si:e  THEN  Out«(3)=Add$ 

2850  OUTPUT  SPathjOutXl) ,0ut$(2),0ut«(3) 

2840  NEH  I 

2870  IF  Iend=Size  THEN  2970 

1880  OUTPUT  sPath; Inbl. lnb«, Inb«, Inb», InbJ, InbJ, Inb«, Inbl, lnb« 

2890  Npg=NpgM 

2900  Tll=Add$iTitIe«!.He3«VALI(Npn) 

2910  htart=lend»46 

2920  Iend=l5tartt44 

2930  IF  Iend>Size  THEN  Iend=SizE 

2940  OUTPUT  »Path;Tl$,Inb»,St$,Inb«,Hel» 

2950  OUTPUT  fPath;He2»,lnbJ,St«,Inb«,Inb$ 

2960  GOTO  2760 

2970  ASSIGN  fPath  TO  I 

2980  HASS  STORAGE  IS  "!HP82901. 700,0" 

2990  BEEP  1302. 08,. 3 
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INPUT  "Do  you  Dish  to  have  a  PLOI  of  the  nunbers  ?  iy/NI',(lns« 

IF  SnsiO'V'  THEN  3200 

PRINT  CHR«(12) 

PRINT  Title*, ■  ","  ■/SOURCE  DflTfl  PRE  :  'jNI 

PRINT 

till     PLEASE  INSERT  fl  DISC  INTO  DISC  DRIVE  1' 


PRINT 
PRINT 
PRINT 
PRINT 


n  BDflT  FILE  KILL  BE  CREATED  ON  THE  DISC.  THIS  FILE  CAN  THEN" 
'BE  PLOTTED  USING  THE  GENERAL  PLOTTING  PROGRAN" 
INPUT  "PLEASE  ENTER  A  UNIQUE  NAHE  FOR  THE  BOAT  FILE  ",Bif« 
MASS  STORAGE  IS  ■:HPB290 1,700,1' 
CREATE  BDAT  6d«, (2l3i2el,9 
ASSIGN  sPath  10  8d« 
Foriat$="(1Z.DDDE,»" 
OUTPUT  fPath  USING  Foriat«i6(ll 
ASSIGN  fPath  TO  I 
iiaSS  STORAGE  IS  ";HP82901, 700,0" 
PRINT 

PRINT  "BDAT  Tile;  "iBd«,Si:e;"  points", "Data  foriat:  BZ.DDDE" 
DEALLOCATE  Gil) 
BEEP  5126. 94,. 1 
PRINT  "     END" 
SUBEND 


SUB  Dec(I,Z(»),Z$(l),Size) 
OPTION  BASE  I 

THIS  SUBROUTINE  ENSURES  THAT  ALL  VALUES  HAVE  THE 
CORRECT  FIELD  LENGTH  FOR  OUTPUT 

Ii=I 

A=0 

Dec«="0" 

Per«="." 

Z«(31="" 

Z«(4)="" 

LDng=5 

FOR  J=l  TO  2 

Z2=ZIIi,Ji 

IF  Z:<0  THEN  Z2«=VAL$(DR0UND(FNRoundlZ2,3l,5ll 

IF  Zz<0  THEN  3550 

Z!«=VAL»iDR0UND(Z2,5ll 

IF  Zz<l.E-5  THEN  Z:«=". 0000000" 

IF  Zz<1.0  THEN  3600 

GOTO  3650 

Z»(JtA)=Z:$ 

Lonq=6 
NEXT  J 

IF  liOI  OR  IM5>Size  THEN  3740 
A=2 

Ii=IM5 
GOTO  3390 

FOR  ZZ<0 

IF  LEN(Zz«)>Long  THEN  Zz«=Z2(tl,LongI 

IF  LEN(Zz$)<Long  THEN  Zz»=ZzHDec« 

IF  LENIZzllOLom  THEN  3570 
SOTO  3480 

FOR  0(ZZ<1 

IF  LEN(Zz»)>=Long  THEN  Zz«=Zz«[l,Long-lI 

IF  LENtZzlKLong  THEN  Iz*=Dec»iZz» 

IF  LEN(Zz$K!Long  THEN  3710 
GOTO  3480 

FOR  ZZ/1 

P«="N" 

FOR  F=LEN(Zz«)  TO  1  STEP  -1 
IF  ZzltF,Fl="."  THEN  P*='y 

NEXT  F 

IF  P«="N"  AND  LENIZz«)<Long  THEN  Zz«=Zz«iPer$ 

IF  Zz<l.Et5  AND  LEN(Zz«)<LonQ  THEN  Zz«=Zz»!.Dec« 

IF  Zzn.Et5  AND  LEN(Zz»(>Long  THEN  3710 
GOTO  3480 
SUBEND 
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3750 

3760 
377(1 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
3990 

4000 

40io 
4020 
4030 
4040 
4059 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
43B0 


DEF  FNRound(Z,Hpl 
Nu=INrailO'Npt.5) 
Nu=Nu/(IO'Np) 
RETURN  Nu 

FNEND 


SU8  6ra(  set(P2,l5cale) 

GINIT  " 

GRAPHICS  OFF 

VIEKPORT  0,130,15,100 

PEN  1 

LDIR  PI/2 

LORG  6 

N8VE  0,57 

LABEL  'pressure  (psiq)" 

LDIR  0 

L0R6  4 

HOVE  65,15 

LABEL  'IIliE  (s)" 

CSHE  3. 7, .5 

MOVE  7,15 

LABEL  -O' 

(IQVE  128,15 

LABEL  Tscale 

LORG  2 

HOVE  0,20 

LABEL  '-5' 

nOVE  0,98 

LABEL  P2 

VIEKPORT  7,130,20,100 

KINDOH  0,T5cale,-5,P2 

FRAME 

A(ES  1,1,0,0,5,5,2 
SU8END 


SUB  Scr  plQtieit),Nuidata) 

BPTIOR  BASE  1 

Er«="N" 

LINE  TYPE  1 

PEN  1 

ALPHA  OFF 

GRAPHICS  ON 

HOVE  S(1.1),6U,21 

FOR  1=1  10  Nmdata 

IF  GII,1)=0  AND  I>1  THEN  4280 

DRAN  Gil, 1), Gil, 2) 

NEXT  I 

NOVE  0,0 

IF  Er«=-Y"  THEN  4360 

NAIT  5 

INPUT  ■     ERASE  THIS  PLOT  (Y/NI?",Er$ 

IF  Er«<)'Y"  THEN  4360 

PEN  -1 

GOTO  4240 

GRAPHICS  OFF 

ALPHA  ON 
SUBEND 


ROUNDING  FOR  7(0 


1X3 

APPENDIX  C 
System  Performance 

Control  System  Performance 

The  ability  of  the  control  system  to  provide  a  constant  power 
delivery  (superficial  heat  flux)  to  the  test  element  is  illustrated  in 
Fig.  CI.   It  will  be  noted  that,  in  spite  of  the  substantial  change  in 
wire  temperature  (and  the  rapidity  at  which  this  change  is  possible), 
with  associated  change  in  wire  resistance,  the  control  system  provides 
constant  power  delivery. 

Pressure  Measurement  System  Performance 

The  ability  of  the  pressure  measurement  system  to  yield  an  accurate 
pressure  history  of  an  experiment  depends  on  proper  setup  for  test 
conditions.   Best  accuracy  requires  flush  mounting  of  the  pressure 
transducer  to  the  test  section.   However,  this  was  not  possible  due  to 
the  necessity  of  providing  cooling  capability  for  the  transducer  with 
the  cooling  adaptor.   The  choice  of  the  amplifier  time  constant  also 
affects  the  accuracy,  and  a  dramatic  effect  on  the  voltage  signal 
measured  for  a  given  pressure  drop.   Use  of  the  short  and  medium  time 
constant  settings  resulted  in  signal  decay  during  the  pressure 
transient,  while  use  of  the  log  time  constant  introduced  the 
complication  of  signal  drift  (Fig.  C2) . 

The  effect  of  temperature  on  the  measure  voltage  signal  is 
illustrated  in  Fig.  C3.   The  figure  compares  measurements  at  an  ambient 
temperature  T^  =  27°C  to  measurements  at  T  =  100°C.   For  the  long  time 
constant  setting,  the  drift  is  adversely  affected  by  Increasing 
temperature  (no  significant  difference  in  behavior  was  found  in  the 
medium  and  short  time  constant  cases) .   Additionally,  as  seen  in  Fig. 
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C4,  there  was  little  consistency  in  the  value  of  the  drift  beyond  the 
individual  run.   Therefore,  the  analysis  program  determined  the  average 
drift  for  each  run  and  corrected  the  data  for  this  drift,  then 
translated  the  data  from  voltage  units  to  pressure  units  based  on  the 
observed  pressure  drop.   Using  the  same  data  as  Fig.  C3,  these  two  steps 
are  illustrated  in  Figs.  C5  and  C6,  comparing  the  results  from  a  test 
with  significant  drift  to  the  results  from  a  test  with  zero  drift. 
Figure  C7  also  illustrates  the  effect  of  ambient  temperature  on  the 
pressure  transient,  comparing  results  for  T  =  95°C  to  results  for  T  = 
lOO-C. 
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FIG.  C7.   Effect  of  ambient  terapci-acure  on  a  0.377  to  0.101  HPa 

pressure  drop.   The  dashed  Hue  is  the  pressure  history 

for  T   =  95''C.   The  solid  line  is  the  pressure  history 

a 
for  T  =  100°C. 
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APPENDIX  D 

Selected  Listings  of  Experimental  Pressure  Data 
Representative  runs  for  series  1-8  (Table  5,2) 
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APPENDIX  E 

Selected  Listings  of  Experimental  Temperature  Data 

First  run  in  each  of  series  1-8  (Table  5.3).   THETA  is  the  platinum 
heater  temperature  less  the  ambient  temperature.   The  NLISSELT  NUMBER  is 
based  on  the  superficial  heat  flux. 
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ABSTRACT  ■.■':', 

An  experimental  Investigation  of  boiling  initiation  during  pressure 
transients  has  been  made.   A  platinum  wire  heating  element  of  dimensions 
0.25  mm  diameter  and  9.6  cm  length  was  immersed  In  distilled  and 
degassed  water.   Maximum  pressures  from  0.377  to  1.48  MPa  were  applied 
with  the  water  and  test  element  temperatures  at  100°C  before  heating  the 
platinum  wire  to  160°C,   The  wire  experienced  constant  superficial  heat 
fluxes  of  approximately  0.42  MW/m^  while  the  system  was  subjected  to 
near-exponential  decompressions  to  atmospheric  pressure.   Pressure 
reduction  periods  were  approximately  4  and  6.6  s,  and  pressure 
transients  were  reproducible  over  most  of  the  pressure  drop.   Boiling 
initiation  times  provided  conditions  for  boiling  initiation  which 
depended  on  the  pressure-temperature  history  of  the  test  element  and  the 
surrounding  fluid.   These  conditions  tended  to  be  overpredicted  by  a 
model  accounting  for  the  deactivating  effect  of  pre-pressurization  on 
potential  nucleatlon  sites.   Reactivation  of  nucleation  sites  and 
recovery  of  the  test  element  temperature  to  steady-state  were  seen  to  be 
strongly  affected  by  increases  in  the  maximum  pressure  applied. 


